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ABSTRACT 


Using an earlier developed mathematical model applicable for 
laminar flow of an incompressible Newtonian fluid, radially inward 
between parallel co-rotating disks, performance curves for an uncon- 
ventional type of turbomachinery are presented. Parabolic and uniform 
input velocity profiles with full admission of the working fluid to 
the rotor is assumed. Performance data for the general flow problem 
require that three dimensionless input parameters be specified, a 
Reynolds number, a flow rate parameter, and a peripheral tangential 
velocity component parameter. The output of the mathematical model 
is the dimensionless characteristics torque, total pressure, and 
efficiency at specified radial locations. 

Using this output data individual performance characterise cs 
are plotted, functions of the Reynolds number and the flow rate 
parameter. Each is a separate constant value curve for a specified 
exit radial location and peripheral tangential velocity parameter. 
Composite maps are also presented showing the characteristi cs of all 
three parameters on a single graph. 

An example is presented which shows how the performance maps 
are used. Using three different working fluids the results, in dimen- 
sional quantities, indicate the specifications and characteristics of 
a multiple-disk turbine. 



All results in this study are for a mathematical model having 
no losses external to the turbine rotor. The performance maps show, 
however, that the multiple-disk turbine is competitive with conven- 
tional bladed turbines in areas of power generation where low 
horsepower and high speeds are required. 
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T 
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r ft 
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radial component of velocity, (dimensionless) 
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v = — tangential component of velocity, (dimensionless) 
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0 
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8,5 dummy variables 



o 


inner radius of disks 
outer radius of disks 


xi 



Chapter I 


INTRODUCTION 

CONVENTIONAL TURBOMACHINERY 

Conventional turbomachi nery utilizes the energy of the working 
fluid to develop pressure, reaction, or impact forces on a series of 
blades or vanes. The blades are generally airfoil shapes, elaborately 
contoured, and precision machined to specifications determined by 
design operating conditions to induce maximum momentum transfer with 
minimum energy losses. 

Methods of analysis used for aircraft wings can be applied to 
study the fluid flow past conventional turbine blades. The flow be- 
tween adjacent blades is complicated and exact solutions of the 
governing equations for the flow process are complex, and a viscous 
flow analysis is difficult. The viscous effects are, therefore, 
usually ignored in analytical work, reducing the problem to an analysi 
of in viscid flow past an airfoil. Empirical correction factors are 
introduced later to account for viscous effects. 

The losses due to viscous shear forces become severe as the 
diameter of bladed rotors and blade spacing decrease. As the turbine 
size is reduced the volume flow rate of the working fluid is reduced, 
therefore there is a reduction in specific speed and efficiency. As a 



result conventional bladed turbines have low efficiencies in situations 
requiring high energy transfer with small volume flow rates. The pri- 
mary advantage of bladed turbines is a high power output-to-wei ght 
ratio and high efficiencies for large units. 

UNCONVENTIONAL TURBOMACHINERY 

One type of unconventional turbomachinery uses fluid viscous 
effects as a mechanism for the transfer of energy. The rotor consists 
of a series of flat, circular disks mounted parallel on a shaft. The 
disks are attached at small intervals along the shaft, and each disk 
has an opening at an inner radius for the exhaust of the fluid. The 
rotor is enclosed with a housing which contains a stator to provide 
the working fluid with a velocity vector that is nearly tangential 
relative to the outer periphery of disks. As the fluid element moves 
in a spiral path to the exhaust port, the viscous shear effects, in 
the tangential direction, tend to draw the disks along in the same 
direction of the flow. A schematic diagram of the fluid flow is repre- 
sented in Figure 1-1. 

Since the geometry of the disk turbine is simple a number of 
analytical and numerical approaches can be used in the flow analysis. 
Any analysis, however, must include the viscous properties of the fluid 
since this property provides the motive power for the multiple-disk 
turbi ne. 

The disk turbine has several inherent advantages. Construction 
costs are relatively low because disks can be produced in large quanti- 
ties from a stamp press; extensive precision machinery is not required. 




HOUS/NS SUPPORT, ROTOR BEARINGS AND BEARING 
SUPPORTS , AND ROTOR TO HOUSING SEALS OMITTED 
FOR CLARITY. 


Figure 1-1. Schematic Diagram of Turbine Using a Disk Rotor 
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Maintenance can be held to a minimum due to the simplicity of the device 
and the ease of replacing damaged or worn disks. The disk turbine also 
allows high operating temperatures to be maintained without structural 
failure. A wide range of working fluids is also possible. The pri- 
mary advantage of the multiple-disk turbine is its ability to produce 
low horsepower, with high efficiency at high speeds, with a very low 
noise level. A typical disk rotor is presented in Figure 1-2. 

BACKGROUND 

The development of multiple-disk turbomachinery dates from the 
turn of the century to the present. Nikola Tesla was granted the IJ.S. 
Patent covering both multiple-disk turbines and pumps in 1913 [1]. In 
his patent application, Tesla described the motive forces of his 
machinery as being dependent on the fluid properties of viscosity and 
adhesion. Several successful turbines were designed and built by 
Tesla [2,3] but they were considered not commercially feasible at the 
time. Lack of adequate materials, difficulties in handling inherent 
high speeds, and heat dissipation problems were the reasons for multiple 
disk devices not being further developed. 

In the succeeding years a number of experimental and analytical 
investigations confirmed the feasibility of multiple-disk turbomachinery 
[4-11]. Since 1955 qualitative analyses have been published which used 
either hydraulic friction factor, empirical information, or assumed 
velocity distributions to estimate performance [12-16]. Although these 
results are not sufficient for design studies, they confirm the feasi- 
bility of multiple-disk devices. 




In 1962 Breiter and Pohlhausen [17] presented rigorous approach 
to laminar flow of a viscous, incompressible fluid between two parallel 
rotating disks operating as a multiple-disk pump. By assuming a very 
small through-flow rate they linearized the governing Navier-Stokes 
equations and obtained a similarity solution. Near the inlet the 
linearized equations did not apply, and it was necessary to solve the 
non-linear equations numerically. The solution of the linearized 
governing equations, valid for small flow rates and radial stations 
well removed from the entrance region, describes the asymptotic flow 
field. The solution of the non-linear equations represents realistic 
inlet conditions and is useful for the radially outward flow case. 

An investigation by Matsch [18] in 1967 considered radially 
inward flow between rotating disks, which models the multiple-disk 
turbine. Reduced Navier-Stokes equations were solved, for both partial 
and full admission of fluid, by an iterative method. The first approx- 
imation was the solution to the linearized form of the governing 
equations. The second order approximation was obtained from the dis- 
carded non-linear terms using the results off the first approximation. 
The procedure was repeated until convergence was obtained, and a 
functional solution through the third approximation was found. The 
solution only specified the dimensionless volume flow rate and Reynolds 
number as parameters, resulting in an insufficient number of parameters 
compared to the number needed for a general solution. The results are 
only valid for the asymptotic flow region occurring at a radial station 
interior to the inlet region. The asymptotic solution does not provide 
sufficient information for the calculation of actual turbine perform- 


ance. 
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A sufficiently complete problem statement was formulated from 
the Navier-Stokes equations by Boyd [19,20] for radially inward flow 
between rotating disks. Boyd introduced a third parameter, a peri- 
pheral tangential velocity component, which may be greater or less than 
the disk peripheral velocity. The reduced Navier-Stokes equations were 
solved for full admission of the working fluid, and an inverse finite- 
difference technique was used to obtain a numerical solution to the 
problem statement. The solution becomes the asymptotic solution at 
interior radii, and is in full agreement with Matsch [18]. The solu- 
tion was the first of sufficient generality to allow the analysis of 
turbine performance over a wide range of flow parameters. 

Adams [21] conducted an experimental investigation which was 
useful in the testing of theory. The test fixture consisted of two 
narrowly spaced disks which were attached to a shaft. The housing 
included supply nozzles which approximated the full admission of fluid 
at the outer periphery. The static pressure distribution was measured 
for the following ranges of the flow governing parameters: 

1.75 < N re < 8.25 

- 0.50 < U < - 0.085 
o 

0.25 < V < 1.75 

o 

where N RE is a Reynolds number and 

rl/2 

U Q = 2 u(l ,z) dz (1-1 ) 


0 
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rV2 

V Q = 2 v ( 1 ,z) dz . (1-2) 

» 

0 

U and V are the mean values of the radial and tangential velocities 
o o 

at the outer periphery of the disks. 

The experimental and calculated pressure distributions compared 
well with the greatest divergence noted for > 6.0. Boyack [22] 
indicates that divergence of results at high Reynolds numbers may be 
due to order-of-magni tude reduced governing equations becoming less 
representative of the actual flow. Other investigators [15,19,20] have 
suggested that the onset of turbulent flow exists in the region. How- 
ever no evidence of turbulent flow was found by Adams for the ranges 
of U , V , and N nlr investigated. 

In 1969 Boyack [22] developed an integral method for the three- 
dimensional, non-boundary-layer flow, which occurs for laminar 
incompressible Newtonian fluid between co-rotating disks. Boyack's 
results agree very closely with the results obtained by Boyd's finite- 
difference method [19] and Adams' experiment [21]. 



Chapter II 


MULTIPLE-DISK TURBINE ANALYSIS 
PRESENT STATUS OF ANALYSIS 

The general laminar flow of an incompressible, viscous fluid 
between co-rotating disks is described by the Navier-Stokes equations, 
subject to the no-slip boundary conditions at the disk surfaces. The 
fluid may be supplied at the outer periphery in many ways. Matsch and 
Rice [23,24] obtained solutions for the limiting cases of potentia' 
flow and creeping flow between co-rotating disks with partial admis- 
sion at their outer periphery. Using order of magnitude arguments to 
simplify the governing equations, Matsch [18] obtained an iterative 
solution for the asymptotic flow field, with full admission. 

The asymptotic solution, however, is not sufficient for the 
calculation of actual turbine performance, because the asymptotic flow 
usually occurs radially inward and away from the entrance. Boyd [19] 
formulated a sufficiently complete problem statement which included 
the entrance region as well as the asymptotic flow. The numerical 
method of solution as presented by Boyd is sufficiently accurate to 
calculate parameters such as torque, power, and efficiency for the 
rotors of multiple-disk turbines. 

Boyd's numerical procedure [19], although accurate, is not 
efficient with regard to time needed on a digital computer to calculate 
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numerous flow cases for the plotting of performance maps. Boyack [22] 
developed an integral solution for the same problem presented in refer- 
ence [19], which is accurate and efficient. The computation time, on 
the digital computer, of the Boyack integral method is one-fiftieth 
(1/50) of that required by the Boyd numerical procedure. 

SUMMARY OF BOYACK 1 S ANALYSIS [22] 

The solution presented by Boyack is necessary for the mapping 
of disk rotor performance for full admission, incompressible, laminar 
flow fields. The method is a forward-stepping procedure which satis- 
fies the integrals of the governing differential equations, boundary 
conditions, and governing differential equations at every radius. 

The Navier-Stokes equations are sufficiently general to 
describe the laminar flow of a Newtonian fluid. With the continuity 
equation the result is four general flow equations. The following 
assumptions are made: 

1. constant viscosity of the working fluid, 

2. constant density of the working fluid, 

3. steady flow, 

4. negligible body forces, 

5. full admission of the working fluid at the outer periphery 
of the disks. 

Using the coordinate system shown in Figure 1 1 - 1 , and de-dimensionalizing 
the momentum and continuity equations they become 



0 


Disks attached at 
2 = + h/2 to shaft 
and rotating with 
angular velocity a. 
Shown is disk at + h/2. 



Figure II-l. Definition of Reference Coordinate System 
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where 


N = P h ^ 
RE 

p 


Using the order-of-magnitude arguments presented by 
h << r , equations ( II-l ) , (I I -2 ) , ( 1 1 - 3) , and (II -4) reduce 


u9u + w9u vf_ _ 1 9 2 u 

9r 9z r 9r N R£ 9z 2 

u9v + w9v + uv _ 1 9 2 v 
9r 9z r N R£ 9z 2 

^0 

9z 



3w 


= 0 


9z 


(II-D 
(II -2 ) 

( 1 1-3) 

( 1 1 -4) 

( 1 1 - 5 ) 

Matsch [18], 
to 

(1 1-6) 

( 1 1 -7 ) 

( 1 1 - 8 ) 


(II -9) 
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The pressure can be assumed to be a function of the radial coordinate 
only. 

If entrance, exhaust, and boundary conditions are assumed 
symmetric about the r - 0 plane midway between the disks, they are 

u ( 1 ,z ) = u o (z) 

v(l,z) = v 0 (z) 0 < z < 1/2 v 1 1 - 1 0 ) 

p(r i ,z) = 0 


r ,0) = 0 u(r,l/2) = 0 

9z 

^{r,0) = 0 v(r,l/2) = r r. < r < 1 (11-11) 

3z i - - 

w(r,0) = 0 w(r,l/2) = 0 

Equations (II-6), (II-7), and (II -9 ) are a system of three 
partial differential equations with four unknowns; velocity components 
u, v, w, and pressure gradient ^ . Another equation is needed. 

Boyack [22] obtained this auxiliary equation by integrating equation 
( I I -6 ) three times with respect to z and using the appropriate boundary 
conditions. The result is 
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To produce the integral equation for the solution, equations 
( 1 1 -6 ) , ( 1 1 -7) , and ( 1 1 -9) , are integrated over the interval 
0 < z < 1/2 and yield 


1 

r 

rl/2 

u 2 dz 

r 

r dr 

i 

l 

) J 

0 


dz 


- 1&+ J_3M< r ,l/2) 
2 dr 9z 


(H-13) 


]_ d_ 
r 2 dr 


1/2 


u v dz 


1 



^r,l/2) 

dZ 


rl/2 

(r u) dz 

0 



2r 


(11-14) 


(11-15) 


Equations (11-12), (11-13), (11-14), and (11-15) constitute the inte- 
gral equations Boyack uses in his investigation. 

Nth order polynomial forms were chosen for the radial and 
tangential velocity profiles. 


N+l 

u(r,z) = l (a 
i=l 

N+l 


v( r,z) = l (b 

j=l 

where a^ and bj are functions of r. 
for continuity be 


i. ,) 2 i_1 . (11-16) 

j. ,) Z+ 1 (11-17) 

It follows that a representative 


w(r,z) 


N+l 

y 

i = l 




i 


(11-18) 


where the prime denotes differentiation with respect to r. 
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Equations (11-16), (11-17), and (11-18) are substituted into 
the integral equations (11-12), (11-13), (11-14), and (11-15), to 
obtain 


dR = _ 
dr 


12 U 


'RE 


N+l N+l 
- + 24 J l 
r i=l j=l 




(i + J - 1 ) ( i + j + 1) 2 


(i+j+D 


(H-19) 


where 


Q . . = a . , a a . , 
w i,J i-l J-l i J-l 


a i-l 


+ a i-i 


b. , b. , 

i -.1 


( 11 - 20 ) 


and 


'Y 7 Vl a j-l + a ^-l a j-l * r ^ a i-l a i- 
j=l 


1 ' b i-1 b ,i - 1 


(i + j - 1) 2( 1+j -’> 


i dp i T 1 n-n a i-i n 

_ _ 2 7 -• o'\ o 


2 dr N re i = l 2 


IT=2 T 


(H-21) 


N y ] N y ] a i-l b j-1 + a i -1 b j-l + r (a i-l b j-l ^ 

(i+j-1 ) 


l l 

i' = l j = l (i + j - 1) 2 


i b i-i 

I- 


n r E 1-1 2 


( 11 - 22 ) 


N+l 


I 

i = l 


a i-l 


(i) 2 


UT 


u j> 

2r 


(11-23) 


Matsch and Rice [25] found that highly inflected velocity pro- 
files occur in the solution. To represent these profiles Boyack had to 
use polynomials of fourth and higher orders. The integral equations 
(11-19) through (11-23) and the boundary conditions are insufficient 
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to determine the coefficients and b. , so additional relations were 
needed. The additional relations between the a^ and b^ were obtained 
by requiring the differential forms of the momentum equations, (II -6) 
through (II-9), to be satisfied in the flow, and yield the correct 
number of equations to allow solutions for a. and b.. 

The additional equations are 


from (1 1-6) 


N+l N+l 

l 1 

i=l j=l 


a . a . 
i J 


!i + a: 

r 1 


j a . b . b . 
J . _! i 


i + 1 


,i+j-2 + d£ 
dr 


, N+l . , 

I (i - 1 ) ( i - 2) a. z 1 " = 0 

n re 1=1 


(11-24) 


from (II -7 ) 


N+l N+l 

1 1 

i=l j=l 



fa. 1 

-L+ a! 

a. b'. - 

i J 

tr ’J 

1 

N+l 

l (T ' 

n re 

i = l 


j b . a . b . ] 

~l ^ 1 1 


i + 1 


„ _] 1 


.i+j-2 


i -3 


= 0 


(11-25) 


The coefficient of z J is zero, so for every value of j, the compata- 
bility relations are obtained from the above equations 


b 2 2 a 9 . 

a a 1 ^+ ^ = 

0 0 r N re dr 


j = 0 


(11-26) 
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A A * 

[ a i + a ,' 

U-l) , b i b j-ll 

a i a j-l ' 

- 


(HI) J* 1 r J 


(j+1)(j+2) a. ? 

= 0 j = 1 ,2 , — (II-27) 


a h 1 - 
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1 

JQ 

03 

H 

X 

i 

i D j-1 

lr ’j 

(HI) “J' 1 r J 


(j+1)(j + 2) b. p 

j = 0,1 ,2, — (11-28) 

N RE 


Using the boundary conditions on the radial and tangential 
velocity components at z = 0, Boyack simplified the compatibility equa- 
tions, (11-27) and (11-28) From equations (11-27) the resulting 
equations are 


2 b n b o 

a 9 a' + a^ a~ 2 

2 o o 2 


+ 


12 a, 


RE 


= 0 


(11-29) 


a 4 a i + a 2 a 2 + a o a 4 


2b o b 4 


b 2 2 


2 

3 


I + a 2 

r 


- 4 



a 4 


30 a g 



0 


(11-30) 


and from equations (11-28) the equations 


a b 2 b 9 

a 0 b ° + - 2 ~ A = 0 


N 


RE 


(11-31) 
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a 2 b i + a o b 2 + (a o b 2 


+ a 2 b o> - 2 


+ a i 


12 b . 


'RE 


= 0 (11-32) 


a 4 b o + a 2 b 2 + a o b 4 + 


(a 4 b o + a 2 b 2 + a o b 4 } 


2 

3 





(11-33) 


The equations Boyack used for the eighth order approximation 
to the exact solution are: the integral equations (11-21), (11-22), 

and (11-23), the simplified compatabi 1 i ty equations (11-29) through 
(11-33), and the boundary conditions which were not used in the pre- 
vious derivations, u(r,l/2) = 0, and v(r,l/2) = r. 

At the outer periphery of the disks, boundary conditions on 
a i and b i are required. If the values of a. and b. are specified at 
the inlet, this corresponds to specifying the distribution of u (z) 
and v Q (z). To obtain the initial values a i ( 1 ) and b.(l), Boyack fit 
to the inlet profiles a general polynomial representative for radial 
and tangential velocity components. The number of points which can 
be fit is pre-determined with the selection of either the fourth, 
sixth or eighth order approximation to the exact solution. The result 
is a system of linear, algebraic equations which can be solved for the 
unknown a^ (1) and b . ( 1 ) . 

A numerical method was used by Boyack [22] to solve his simpli- 
fied, initial value problem consisting of a system of ordinary 
differential equations and appropriate initial conditions. For each 
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radial step, Boyack found a solution could be obtained more efficiently 

by a predictor-corrector method. Although this method has very rapid 

computational time and the readily estimated error can be used as a 

basis for choosing step size, the solution must be known at r 

r r ,, and r in order to find the solutions at r , , . Thus the 
n-2 n-1 n n+1 

predi ctor-corrector method is not self-starting. For this reason, 
a Runge-Kutta-Gi 11 method was used to provide the starting solution 
and an Adams-Moulton predictor-corrector method was used following the 
starting solution. These methods are well documented by Boyack [22]. 

Using a small step size the solution was started at r = 1. 

For each step forward the truncation error, E^, was calculated and 
compared with preset allowable error limits, E and E . . If 

I [id A 1 1 1 1 f i 

E, > E the step size was halved and the Runge-Kutta-Gil 1 method was 
K max 

used to re-start the solution. If E„ < E m . , the step size was doubled 

i\ rni n 

and the solution was re-started. If E . < E„ < E_ . the step size 

ill 1 n l\ HldX 

remained constant and the A-M method was used to continue the calcula- 
tion . 

When comparing results for fourth, sixth and eighth order solu- 
tions with each other, Boyack found that the results for the fourth 
order solution were very poor due to the u velocity component becoming 
highly inflected at some radii. The results of the sixth and eighth 
order solution were in good agreement, indicating the solution approaches 
the correct solution as the order increases. Comparing the eighth order 
with the results obtained by Boyd [19] and Adams [21] the agreements 
were excellent. A complete outline of Boyack's digital computer pro- 
gram is given in Appendix A. 
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TURBINE PERFORMANCE PARAMETERS 
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A complete derivation of pertinent turbine performance para- 
meters may be found in reference [19]. By taking a control volume 
consisting of the fluid between the disk inner and outer radii and 
applying a moment of momentum equation, the net torque on the turbine 
shaft is found. If all external rotor losses are neglected, the shaft 
torque is the negative of the control volume torque, and can be written 
as 


T 

s 


• 1/2 


4tt r‘ 


<u. 


v o> 


dz - 4tt r 




1/2 


(Ui 


v i } 


dz 


(11-34) 


The total pressure is obtained from Bernoulli's equation and can be 
written as 


Pt 



1/2 

( U n + V n ) U r > dZ 
0 0 0 


( IT-35) 


The two above parameters can be expressed in terms of a^ and b. if 
equations (11-16) and (11-17) are substituted into equations (11-34) 
and (11-35). As a result of the substitutions the following expres- 
sions are obtained 
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i = l j=l 
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(11-36) 


1_ a 0j ( a 0j a 0k + b 0j b _°k) 
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Pt = Po + 
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The efficiency of the turbine is expressed as the ratio of actual 
shaft work to the ideal shaft work, where the ideal shaft work is 
derived from the energy equation and the mass flow. The efficiency 
can be expressed as 


n 


T 

s 


2 " U o p t 


(11-38) 



Chapter III 


PERFORMANCE MAPS 

The primary objective of this study was to map the performance 
characteristics of the multiple-disk turbine. A designer can utilize 
the data shown in the maps to produce a wide variety of design charts 
in terms of dimensional quantities for a particular design study. 

BOYACK'S COMPUTED RESULTS 

For each combination of input parameters, U , V Q and N RE a 
summary sheet is produced which tabulates the dimensionless turbine 
performance parameters total pressure (PT), pressure (P), torque (T), 
and efficiency (n) at specific radii R. In the summary R is the value 
of dimensionless inner radius of the turbine; therefore, the summary 
sheet covers an infinity of multiple-disk turbines. It was decided 
that the three dimensionless performance parameters, T, PT, and n could 
best describe the turbine characteristics. PT was chosen in lieu of P 
because P represents the pressure drop of the flow field within the disk 
region only. The value PT accounts for a pressure drop across the in- 
let nozzles as well as the disk region. The pressure drop across the 
nozzles does not remain constant; therefore the PT values provide a more 
adequate representati on of the flow field. 
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Three different types of inlet velocity profiles can be assumed 
using the Boyack program, Matsch asymptotic, parabolic or uniform. The 
Matsch asymptotic profiles give results only applicable to the asymp- 
totic flow region. The results obtained using parabolic and uniform 
velocity profiles are adequate to determine turbine characteristics 
because the solution of these profiles consider inlet conditions. It 
was decided to concentrate on the parabolic inlet profiles, because it 
is felt that the velocity profiles of the fluid exiting the stator 
nozzles is sufficiently developed to be parabolic in nature. 

The Boyack program requires approximately 30 seconds to compute 
the performance parameters for a given combination of U Q , V Q and Np^, 
using the CDC 6400 computer; making it extremely efficient as compared 
to the Boyd finite-difference method. A sample of the data summary 
sheets used in this study are presented in Appendix B. 

It must be emphasized that all results presented are for a 
mathematical model having no losses external to the rotor. The model 
assumes perfect nozzles, perfect exhaust diffusion, no bearing or seal 
losses, and no windage losses between end disks and the rotor housing. 

GRAPHICAL METHODS 


The three performance parameters, n» T, and PT, were plotted 
as functions of U Q and Np^ at specific radial locations and constant 
V Q . To do this, a set of summary sheets was obtained and initial 
graphs were plotted, n , T, and PT versus U Q were plotted separately 
with Np£ and V Q constant at specific radial locations. Ten discrete, 
uniform radial locations were chosen, 0.1 to 0.9. The data from the 
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initial graphs was then cross-plotted to obtain constant value 
performance curves. To obtain these curves, constant values of n, T, 
and PT were taken from the initial graphs and plotted separately at 
corresponding values on a U 0 versus N RE graph, with V Q and ft again 
constant. These points were later joined to form smooth and contin- 
uous curves of constant values for each of the nondimensi onal turbine 
parameters. 

In a preliminary investigation of parabolic inlet profiles the 


input values for N RE and U Q were chosen as follows: N RE = 2.0, 4.0, 

6.0, 8.0, and 10.0; 0.02 < U Q < 1.0, using an interval of 0.02 between 
0.02 and 0.3, and an interval of 0.1 from 0.3 to 1.0. The value of 


V q - 1.1 was also used in the preliminary investigation. These points 
were selected to obtain, generally, the operating limits of the multiple- 
disk turbine. This initial investigation did indicate the operating 
limits and also areas which required further investigation. 

It was found that for input values of U Q greater than 0.3, at 
all N re values, efficiency was generally less than 20%, and total pres- 
sure extremely high, making this region of little interest in turbine 
design. The curves on the initial graphs were smooth and predictable, 
so it was decided that the narrow interval between the input values of 
U Q could be relaxed. At lower values of N R ^ it was found that n and PT 
undergo greater changes, and additional input values of N RE had to be 
used to determine inflection points. 

Using the preliminary investigation as a guide, further investi- 
gations of parabolic inlet profiles were conducted for input values of 
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V Q = 0.8, 1.0, and 1.3. An investigation of uniform inlet profiles 

was also conducted for an input value of V =1.1. 

o 

Figures 1 1 1 - 1 through II 1-3 are typical performance maps for 
constant values of dimensionless torque, total pressure, and efficiency. 
In addition composite maps were constructed, Figure 1 1 1 -4 , which show a 
reduced number of curves superimposed on a single map for all three tur- 
bine parameters. These composite maps give a clearer understanding of 
the behavior of the turbine parameters with respect to each other. 
Complete sets of performance maps for parabolic inlet profiles with 
V = 1.0, 1.1, and 1.3 and uniform inlet profiles with V =1.1 are 
presented in Appendix C. 

A large number of one-to-one graphs were plotted to verify the 
characteri sti cs obtained on the performance maps. One set of graphs 
is constant value R. curves plotted for corresponding q, T, or PT and 
V Q with U Q and constant. The other set of one-to-one graphs is 

constant value curves plotted for corresponding n and R. with U Q 

and V Q constant. A selected number of these curves are in Appendix C. 

GENERAL CURVE CHARACTERISTICS 

Looking at the performance maps for a fixed V Q , it is evident 

that x] increases as and U Q decrease. The dimensionless torque T 

increases as R. decreases and it decreases as U decreases. PT the 
i o 

dimensionless total pressure decreases as U decreases and it increases 

o 

as R^. decreases. In the region of the maps where U Q is greater than 
0.10, efficiency is low; total pressure and torque are high. Here the 
performance maps indicate that high torque turbines are possible at 
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the expense of low efficiency and high total pressure. This region of 
operation has questionable utility. 

The area of prime interest lies in the region where U Q is less 
than 0.05. Through the range of N^, the T curves have relatively zero 
slope for any given U . The PT curves in this region do not have 
drastic inflections and their values are small. It follows that at a 
particular value of U , p is a function of PT only since 

0 = (III-D 

Zir U PT 
o 

As previously stated T decreases as U Q decreases so this region is ideal 
for the design of small, low powered turbines, because high efficiency 
and low total pressure requirements are easily obtained. 

If a particular value of torque is required, it is possible to 
choose from the performance maps a design value for the exhaust radius 
that will provide maximum efficiency; or if a specific total pressure 
drop is available the exhaust radius can be picked to provide maximum 
efficiency, and the corresponding values of U and N DC . are fixed. There 
are any number of combinations available. 

The first set of one-to-one plots (o versus V ), verify that 
efficiency increases as decreases to a value of about 0.3; then the 
efficiency starts to decrease. The maximum efficiency occurs within 
one-tenth of V Q = 1.0. The other two parameter graphs indicate the 
values of PT and T decreases as R. increases for all V . 
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The other set of one-to-one plots (n versus R . ) shows that the 
efficiency peaks at equal to about 0.3 and = 4.0 for all U Q 

and V . 

o 

The performance maps for the uniform inlet profiles follow the 
same pattern as the ones for the parabolic inlet profiles. During the 
investigation it was found that the polynomial representation for uni- 
form inlet profiles presented in reference [22] was incorrect. A 
suitable representation , which could be solved by the Boyack integral 
method, was not found. Crawford [29], during his analytical design and 
optimization study of a multiple-disk pump, developed a method using 
the Boyd [19] finite-difference technique to obtain a solution in the 
entrance region. When the velocity profiles can be represented by a 
polynomial, without discontinuities, the Boyack integral method is used 
to continue the solution. The profile match-up occurs before the flow 
reaches R. = 0.98 and generally takes 60 finite-difference steps. 

The area of interest on the performance maps for the uniform 
inlet profiles is the region where the input values of U o are less 
than 0.06. The scale of the maps was changed to show only this region. 
While the patterns of the curves remain the same as the parabolic inlet 
profiles, the values of the parameters change significantly. The effi- 
ciency increases approximately 10%, while the total pressure decreases 
and the torque remains about the same for a given 1) . 

LIMITATIONS 

The investigation of the multiple-disk turbine performance para- 
meters disclosed that the turbine was limited in three areas. The first 
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occurs when the inlet profiles are parabolic and the input parameters 
a re U Q < 0.02 ’ N RE > 4, °’ and all V Q . In this region Crawford dis- 
covered that the velocity profiles collapsed, recovering for some 
cases, but not for others. This explains the discontinuity indicated 
on the performance maps for parabolic inlet profiles with V Q = 1.3. 

As N re becomes larger with U Q < 0.02, the velocity profiles show that 
the flow between the disks becomes reversed and high shear forces pre- 
vent the fluid from continuing inward. 

The other two limitations are generally impractical. First 
is the combination of input parameters, U >1.0, all N DC , and all V . 
The efficiency remains extremely low and the total pressure drop 
required to turn the turbine is extremely high. In some cases the 
efficiency remains negative, a finite distance beyond the entrance 
region, indicating a pumping action. 

The other impractical area is the combination of input para- 
meters, V Q < 1.0, all U Q , and all N R ^. V q is defined as the mean 
value of v( 1 ,z) and 

v(l ,z) = 7 ^— ( II 1-2) 

F n 

0 

If v (z) is less than 1.0, the initial tangential velocity of the 
fluid must be less than the tangential velocity of the outer periphery 
of the disks. The disks must do work on the fluid to bring its velo- 
city up to velocity of the disks, and is so indicated in the data 
summary sheets as a negative efficiency. The velocity of the fluid 
attains the velocity of the disks before the flow reaches R.. = 0.95, 
but not before a finite amount of turbine shaft work is lost. 
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While the operation of the multiple-disk turbine in the regions 
described above seem to have little utility, there may be special 
situations where operating in these areas may be useful. 

DESIGN CHARACTERISTICS 


A turbine designer can approach a design problem in many ways. 
Physical constraints and operating conditions dictate which approach 
to design must be taken. The performance maps presented in Appendix C 
are beneficial in visualizing the dependence of n> T, and PT on , 
U Q , V Q , and R^ . Using the maps, a number of design charts in terms of 
dimensional quantities can be produced for a particular study. The 
problem depends on parameters which are fixed, known and unknown. The 
result can be one of a large number of possible combinations of input 
and output data. 

Presented is an example of an initial design study of a 
multiple-disk turbine rotor. The first step is to dimensional ize the 
performance map parameters. From reference [19] the following equa- 
tions for the dimensionless parameters are obtained. 
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h 2 fl 
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(III-5) 


Solving equations (III -3) and (III -4) for dimensional total 
pressure and torque the following is obtained. 


PT = PT p Q 2 r Q 2 


and 


T = T p r 4 h 
o 


Using equations (III -3) to obtain the angular velocity gives 


PT 


PT p f o 2 


1/2 


( 1 1 1 -6 ) 
(III-7) 


( 1 1 1 -8) 


In reference [22] a form of the integrated continuity equations at the 
outer radius is 


Q = 2tt fi h U r 2 
oo 


( 1 1 1-9) 


where Q is the volume flow rate. With the above equations n can now be 
presented as a function of dimensional parameters. 


T fi 

n = — — 

Q PT 

Solving equations ( 1 1 1 -5) for h, the disk spacing, yields 


(III-10) 


h = 


( N Re / 1/2 


(iii-ii) 


Finally the inlet nozzle angle a is defined as 


a = tan 


-1 % 


( 1 1 1 -1 2 ) 


All turbine parameters are now in terms of dimensional quantities. 
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A digital computer program was developed to use the input and 
output results of the modified Boyack program with the dimensional 
total pressure (head) and outer rotor radius available as additional 
input data. The program outputs are the dimensional quantities needed 
for initial design studies. In addition the program computes off design 
point conditions caused by a varying head. A listing of the design 
program is presented in Appendix D. 

In an example presented herein it is assumed that efficiency 
must be at least 80%, the maximum head available is 220 feet, and the 
velocity inlet profiles are parabolic. From the performance maps ' ! t 
is seen that efficiency approaches 80% in the region where = 0.3, 

U Q = 0.02, V Q = 1.1, and = 4.0. The data summary sheets indicate 
that with those input parameters the output characteri sties are 
n = 79.5, PT = 1.525, and T = 0.152. Using the above nondimensi onal 
parameters, fixing the disk thickness and rotor length-to-radi us ratio, 
the dimensional characteristics were calculated for a given working 
fluid. 

Using liquid hydrogen, liquid sodium, and glycerin, and assum- 
ing they behave as incompressible Newtonian fluids, the dimensional 
operating characteristics of a multiple-disk rotor were studied. 

Figure II 1-5 shows how the design angular velocity (RPM) of the rotor 
and output torque ( FT -LB) vary as the fluid density varies with 
selected outer radii. As the fluid density and outer radius increase 
the RMP decreases while the torque increases at design conditions. 

High power outputs occur in large rotors using dense fluids. 
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The graphs in Figures 1 1 1 -6 through III -8 illustrate how the 
design angular velocity, torque, and disk spacing (inches) vary as the 
design head varies for a given rotor size. When the design head in- 
creases, the disk spacing becomes narrower, so for a given size rotor 
the number of disks must increase to meet the design characteristics. 
Denser fluids require wider disk spacing, and a lesser number of disks. 
The torque and angular velocity increase as the design head increases, 
and it follows that the output power increases. For a fixed size rotor, 
the increase of the design head produces an increased volume flow. 

Figures 1 1 1 -9 through 1 1 1 - 1 1 show how design angular velocity, 
dis< spacing, and torque vary as the design rotor size varies for a 
given head. Here again it is evident that larger rotors require wider 
disk spacing and produce higher torque. From equations ( 1 1 1 -9 ) Q is a 
function of ft, h, and r Q 2 . Therefore, in large rotors volume flow 
becomes significant. 

The figures described above were for design conditions for the 
various parameters mentioned. Figures 111-12 through 1 1 1 -1 4 illustrate 
a typical off-design-point operation. For a given design the fluid de- 
sign head was allowed to vary ±20%. The efficiency of the off design 
operation was calculated using equation (III-10), and did not vary. 

The torque and the power varied ±20% and ±30% respecti vely. The volume 
flow, however, only varied ±10%. The above percent changes are inde- 
pendent of rotor size and working fluid. 

The study presented is an example of how a design problem can 
be approached. For a specific turbine requirement there are numerous 
combinations of parameters which can be varied to obtain an optimum 
design. 








Figure 1 1 1 - 7 . Plot of RPM, Disk Spacing, and Torque 
Versus Head for Liquid Sodium and Outer Rotor Radius 
= 4.0 Inches 



Figure III-8. Plot of RPM, Disk Spacing, and Torque 
Versus Head for Glycerine and Outer Rotor Radius = 
4.0 Inches 






Figure 1 1 1 -9 . Plot of RPM, Disk S 
Versus Outer Rotor Radius for Liqu 
Head = 200 Feet 







Figure III-10. Plot of RPM, Disk 
Versus Outer Rotor Radius for Li q 
Head = 200 Feet 












Figure III-13. Illustration of Off-Design-Point 
Operation for Liquid Sodium. Design R = 4.0 
Inches and Head = 200 Feet 0 









Chapter IV 


CONCLUSION 

SUMMARY 


The results of the Boyack analysis were used to plot performance 
maps for the multiple-disk turbine. The Boyack solution is applicable 
for the conditions of a three-dimensional, laminar flow, of an incom- 
pressible Newtonian fluid, in radially inward flow between parallel, 
co-rotating disks. Parabolic and uniform input velocity profiles, with 
full admission of the working fluid to the rotor, were assumed. Per- 
formance data for the general flow problem required the three 
dimensionless input parameters be specified, N^, U Q and V Q . In this 

investigation N DC , U , and V were allowed to vary over the ranges of 
Kb 0 0 

0.5 < N dc < 10.0, 0.01 < U < 1.0, and 0.8 < V < 1.3. The performance 
parameters, n» T, and PT, at specified exit radii were plotted indi- 
vidually with respect to Reynolds number, N^,and radial input velocity. 


The performance maps in Appendix C enable a turbine designer to 
determine maximum limits of efficiency and performance for a multiple- 
disk turbine. Having used the performance curves to determine a 
particular area of interest, a designer could then refer to the output 
sheets presented in Appendix B to obtain more precise data for calcu- 
lating the turbine operating specifications. The dimensional 
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characteristics such as disk spacing, exhaust radius, and nozzle angle 
can then be determined by using a design program similar to the one 
presented in Appendix D. 


RECOMMENDATIONS 


The constant characteristic values on the performance maps were 
hand plotted using the output data from Boyack's analysis. The process 
was very time consuming. It is recommended that a digital computer 
program be constructed, using the output data on the summary sheets in 
punched card form, to make initial plots of 1J versus n, T, and PT 
with N Dr , V , and R. constant. A "parabolic fit" technique can be used 
to make the plots, since the initial curves are smooth and without 
inflections. The results could be graphical or tabular depending on 
the peripheral equipment available. 

A more sophisticated program could be developed to do the cross 
plotting by interpolating a series of initial output curves. The out- 
put could be graphical or tabular again depending on the peripheral 
equipment available. This approach would result in greater accuracy 
and a more efficient method of plotting performance maps. 


FURTHER STUDY 


Performance curves that are produced for conventional turbines, 
as in references [27] and [28], are expressed in terms of specific 
diameter versus specific speed, where dimensional parameters such as 
pressure head, volume flow rate, speed, and diameter are combined to 
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form the specific parameters. Further investigation into the derivation 
of similar expressions for the multiple-disk turbine would be benefi- 
cial. The output data shown in the performance maps in Appendix C 
could be represented by a single set of performance curves rather than 
three individual sets. 

An additional area for investigation is the region where 
U Q > 1.0, and all V Q and N^. With these input parameters the output 
data indicates that there is pumping action in the outer periphery of 
the rotor. The possibility of a radial inflow pump exists in this 
area of operation. 



REFERENCES 


1. Tesla, Nikola. "Turbine," United States Patent No. 1061206, 

May 6, 1913. 

2. "The Tesla Turbine," World's Work , pp. 543-548, March 1912. 

3. "The Tesla Steam Turbine," Scientific American, pp. 296-297, 

September 30, 1911. 

4. Schroeder, Henry B. "An Investigation of Viscosity Force in Air 

by Means of a Viscosity Turbine," Unpublished thesis, 

Rensselear Polytechnic Institute, New York, 1950. 

5. Leaman , A. B.. "The Design, Construction and Investigation of a 

Tesla Turbine," Unpublished Masters thesis. University of 
Maryland, 1950. 

6. Armstrong, J. H. "An Investigation of the Performance of a 

Modified Tesla Turbine," Unpublished Masters thesis, Georgia 
Institute of Technology, 1952. 

7. Young, S. K. "The Investigation and Analysis of the Tesla Turbine," 

Unpublished Masters thesis. University of Illinois, 1957. 

8. Gruber, E. L. "An Investigation of a Turbine with a Multiple Disk 

Rotor," Unpublished Masters thesis, Arizona State University, 
1960. 

9. Smith, F. J. "An Investigation of a Disk Type Pump," Unpublished 

Masters thesis, Arizona State University, 1960. 

10. Gordon, William. "An Investigation of a Disk Type Compressor," 

Unpublished Masters thesis, Arizona State University, 1962. 

11. Hasinger, S. H. and L. G. Kehrt. "Investigation of A Shear Force 
Pump," ASME Transactions, Journal of Engineering fo r Power, 

Volume 85, Series A, No. 3, July 1963. 

Rice, Warren. "An Analytical and Experimental Investigation of 
Multiple-Disk Pumps and Compressors ," ASME Transactions, Journal 
of Engineering for Power , Volume 85, Series A, No. 3, July 1963. 


12 . 



51 


13. Rice, Warren. "An Analytical and Experimental Investigation of 

Multiple-Disk Turbines," ASME Transactions. Journal of Engi- 
neering for Power , Volume 87, Series A, No. 1, January 1965. 

14. Beans, W. E. "Investigation into the Performance Characteristics 

of a Friction Turbine," Journal of Spacecraft , Volume 3, No. 1, 
January 1966. 

15. Nedl , D. "Eine Theoretische Betrachtung der Reibungsturbomashinen 

von Nikola Tesla," Doktor-Ingenieurs dissertation, Lehrstuhl 
fur Technische Mechanik, Technische Hochschule Aachen, July 
1966. 

16. Vanerous, T. "Rotierende Scheiben for Luftvorwarmer mit 

Geblaseworking," Allegemeine Waermetechnik , 1955. 

17. Breiter, M. C. and K. Pohlhausen. "Laminar Flow Between Two 

Parallel Rotating Disks," ARL, USAF, Dayton, Ohio, March 1962. 

18. Matsch, L. A. "Inward Flow of an Incompressible Fluid Between 

Rotating Parallel Disks," Unpublished Doctoral dissertation, 
Arizona State University, May 1967. 

19. Boyd, K. E. "An Analysis of the Laminar Flow of a Viscous, 

Incompressible Fluid Through a Full Admission Multiple-Disk 
Turbine," Unpublished Doctoral dissertation, Arizona State 
University, 1967. 

20. Boyd, K. E. and W. Rice. "Laminar Inward Flow Between Rotating 

Disks with Full Peripheral Admission," ASME Transactions , 

Journal of Applied Mechanics , Volume 35, Series E, No. 2, 

June 1968. 

21. Adams, R. G. "An Experimental Investigation of the Inward Flow 

of an Incompressible, Viscous Fluid Between Parallel, 

Co-Rotating Disks," Unpublished Doctoral dissertation, 

Arizona State University, May 1969. 

22. Boyack, Brent E. "An Integral Solution for the Three Dimensional, 

Laminar Flow of a Viscous, Incompressible Fluid Between 
Parallel, Co-Rotating Disks," Unpublished Doctoral dissertation, 
Arizona State University, June 1969. 

23. Matsch, L. A. and W. Rice. "Potential Flow Between Two Parallel 

Circular Disks with Partial Admission," ASME T ransactions , 
Journal of Applied Mechanics , Volume 34, Series E, No". TV 
March 1967, pp. 239-240. 



52 


24. Matsch, L. A. and W. Rice. "Flow at Low Reynolds Number with 

Partial Admission Between Rotating Disks," ASME Transactions, 
Journal of Applied Mechanics , Volume 34, Series E, No. 3, 
September 1967, pp. 768-770. 

25. Matsch, L. A. and W. Rice. "An Asymptotic Solution for Laminar 

Flow of an Incompressible Fluid Between Rotating Disks," 

ASME Transactions, Journal of Applied Mechanics , Volume 35, 
Series E, No. 1, March 1968, pp. 155-159. 

26. Dant, Robert E. "Calculated Performance of a Multiple-Disk 

Turbine Operating with Incompressible Fluid," Unpublished 
Masters thesis, Arizona State University, June 1969. 

27. Balj£, 0. E. and D. H. Silvern. "A Study of High Energy Level, 

Low Power Output Turbines," AMF/TD No. 1196, Sund Strand Turbo, 
Division of Sund Strand Machine Tool Co., April 1958. 

28. Balj£, 0. E. and R. L. Binsley. "Performance Prediction: 

Optimization Using Fluid Dynamic Criteria," Norr 4507(00), 
Rocket Dyne, Division of North American Aviation, Inc., 

December 1966. 

29. Crawford, Michael E. "A Composite Solution Method for Analytical 

Design and Optimization Studies of a Multiple-Disk Pump," 
Unpublished Masters thesis, Arizona State University, February 
1972. 



APPENDIX A 


MODIFIED BOYACK PROGRAM 
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PROGRAM PUMP C I NPUT , OUTPUT ) 

PROGRAM PUMP, M CRAWFORD, ARIZONA STATE UNIVERSITY, 1971 
PROGRAM UPDATED FOR CDC 6400, SCOPE 3.3 

THIS IS A MODIFICATION OF THE TURBINE PROGRAM, YIELDING 
MAGNITUDE (NON-DIMENSIONALIZED BY OUTER RAOIUS * OMEGA) 

ANO DIRECTION (DEGREES FROM A RADIAL LINE) OF VECTOR AVERAGE 
VELOCITY AT SELECTED RADIAL STATIONS, BY MEANS OF SUBROUTINE 
VECTOR. 


THIS PROGRAM TREATS THE LAMINAR FLOW OF AN INCOMPRESSIBLE, 
NEWTONIAN FLUID, WITH RADIALLY INWARD OR OUTWARD THROUGH- 
FLOW, BETWEEN PARALLEL, CO-ROTATING DISKS. 

THE SOLUTION METHOO CONTAINS A FINITE-DIFFERENCE SOLUTION SEG- 
MENT AND AN INTEGRAL SOLUTION SEGMENT. 

INPUT VELOCITY PROFILE SHAPES MAY BE PARABOLIC , UNIFORM ,MATSCH 
ASYMPTOTIC, OR ARBITRARY SPECIFIED. 

THE F-D SEGMENT IS USED AS A START-UP FOR UNIFORM INPUT PROFILES, 
WITH MATCHUP TO INTEGRAL METHOD UPON SATISFACTION OF MATCHUP 
CRITERIA. 

POLYS OF 4,6,8 ORDER MAY BE SPECIFIED FOR INTEGRAL SEGMENT FOR 
HIGHER OROER PROFILES, REPLACE SUBROUTINE OERV. THE INTEGRAL 
SEGMENT ALLOWS ORDER REDUCTION TO 4 OR 6 DURING INTEGRATION. 



PARAMETER 

TYPE 

DESCRIPTION 

COLUMNS 


* *** ***** 

* * * * 


******* 

CARO 1 

****** 


LIMP 

INT 

NUMBER INFO STORAGE POINTS 

1-2 

CARD 2 

** **** 


XPNT (I) 

FP 

STORAGE POINTS 

4F15. 0 

CARO 3 

(UNIFORM PROFILE 

CONTROL) 


****** 


KK1 

INT 

NO. F-0 STEPS BEFORE MATCHUP TRY 

1-5 


KK2 

INT 

INTERMED. PRINT INTERVAL , WI THIN 

KK1 6-10 


KK3 

INT 

KK3=0 NO PRINT AT KK2 INTERVAL 
KK3= 1 PRINT 

11-15 


PCTU 

FP 

HA X PERCT DEV, INTEGRATED F-0 PRO- 16-30 
FILE TD POLY CURVE-FIT (U VEL) 


PC TV 

FP 

MAX PERCT DEV, V VEL F-0 POINT TO 
POLY CURVE-FIT POINT 

31-45 


NOTE -- IF COMPARI SISON TESTS FAIL, F-0 SEGMENT CONTINUES 
WITH MATCHUP ATTEMPTS AT 10 F-0 INTERVALS 


C CARO 4 

c ****** 

C NOPNT INT NOPNT=0 PRINT EVERY XPNT INTERVAL 1-5 

r - NOPNT = l PRINT ONLY SUMMARY SHEET 



f\) H* <X> C® N ^ ^ (\j H> 


55 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 

r 



MODE 

INT 

M00E=4 6TH ORDER POLY 

6-10 


INPRO 

I NT 

M0DE=6 6TH ORDER POLY 
M0DE=8 8TH ORDER POLY 
INPRO-1 PARABOLIC INPUT PROFILE 

11-15 


KCMODE 

INT 

INPRO=2 UNIFORM INPUT PROFILE 
INPR0= 3 MATSCH ASYMPT PROFILE 
INPR0=4 ARBITRARY PROFILE 
KCMOOE=0 NO PROFILE ORDER REDUCE 

16-20 


XFIT 

FP 

KCM00E=4 4TH OROER RECUCE 
KCMO0E=6 6TH ORDER REDUCTION 
INTERVAL FOR REDUCTION IN POLY 

21-35 

CARO 5 

****** 


UO,VO,PE 

FP 

PARAMETERS FOR INPR0=1,2 

3F15. 0 

OR 


UO,PE 

FP 

PARAMETERS FOR INPR0=3 

2F15 . 0 

OR 


A UP( I , J) 

FP 

Z COORD, U COMP, V COMP ARRAY, 

3F15. 0 

NOTE 

-- DATA CARDS 4,5 

(MODE/2+1) PTS REQD 

REPEAT AS NEW OATA. USE ONE BLANK 

* * 

CARO TO 

* * 

TERMINATE PROGRAM 

******** 


DIMENSIONS FOR 12TH 

ORDER 

VELOCITY PROFILES AND 31 PRINT STATIONS 


WITH PROFILE VALUES AT 21 VALUES OF 0.025 SPACING 

TO COMPUTE ABOVE 8 TH ORDER, SUBROUTINE DERV MUST BE REPLACEO 
WITH THE ARBITRARY ORDER PROFILE DERV SUBROUTINE 

COMMON QS( 13, 13) ,F (13,13) ,G( 13,13) ,H( 13,13) , AC (13) ,8C(13) , 
Y(12),Y1(12),Y2(12),Y3(12),Y4(12),Z(12),Z1(12), 

Z2(12) ,Z3(12) , 74(12) , YI(12) ,E(12) , XI , XC , X , XI , X 2 , 

X3, X4, OXI, OX,OXM, EMIN, EMA X , KONT , LC , MP, MO , Ml , M2 , 

M3 , H4, MT,NT,RF,UO,VO,OPDR,P,P1,P2,P3,P4,PG1,PG2, 
MOOE,NOPNT , INPRO ,KKONT ,MOIOG,KF ,KFF ,KPNT,SIG, ALPRM, 

ATES (6 4 , 65 ) ,XTES<64) , 

RSV (31) , ETASV ( 31 ) ,PTSV(31) ,PSV(31) ,TQUESV(31) , 

OPORSV (31) , ACO (13) ,ZINF(13) ,ZVUV(21,13> ,ZVW(21,13) , 

INUH ( 3 1 ) ,R.INF (11) , A UP (7,3) , UW ( 21) , V W (21 ) ,RI1, 

Z W ( 21) , WW< 21) , Q (13,13) ,XPNT(31) ,RI2 , IS, LU, NTP, PI ,NTH, 

VHA G (3 1 ) , VANG( 31) 

DIMENSION IFIT<7) 

EQUIVALENCE (X,R) 

*1* PE AO NUMPEP OF PRINT ANO/OP DATA STORAGE STATIONS 
READ 7000, LIMP 
* 2 * READ STORAGE STATION ARRAY 

READ 7005, (XPNT(I), 1=1, LIMP) 

READ UNIFORM PROFILE CONTROL DATA 


READ 7015, KK1,KK2,KK3*PCTU,PCTV 

C *4* REAO °RINT CON TRO L , PROF I LE ORDER , PROF ILE SH APE , RE DUCT I ON ORDER, RED. PT 



ooo 


1 REAO 101, NOPNT,MOOE, INPRO , KC MOOE , XF IT 
IF (MODE. EG. 0) CALL EXIT 

GO TO (205,205,210,215), INPRO 

PARA 8 QL IC INLET PROFILES 
OR 

UNIFORM INLET PROFILES 
READ 102, UO , VO, RE 
GO TO 220 

r * 5 * REAO MATSCH ASYMPTOTIC PROFILES 
210 REAO 103, UO, RE 
GO TO 220 

C 

c *5* REAO ARBITRARY PROFILES AS SETS OF Z COORD, U VEL COMP,V VEL 
C COMP FOR ( MOOE / 2+1 ) POINTS, INCLUDING TWO BOUNDARY POINTS 
c AUP(I,1)=Z COORD AUP ( I , 2) -U COMP AUP(I,3)=V COMP 
215 REAO 104, RE 
LlM=M00E/2+l 

REAO 105, ((AUP(I,J), J=l, 3) , 1=1, LIM) 

220 CONTINUE 
IUNTFM* 0 

IF (INPRO. EQ.2) I UNIFM=1 
C CURVE FIT POINTS 

CALL CFIT ( IF IT ,MOOE) 

0 

C$$s$ CDC 6400 TIME RECORD 
BTIME=SECOND( A) 

C 

C XI INITIAL INTEGRATION POINT 

XI=XPNT (1) 

C FINAL CUT-OFF VALUE INDEPENDENT VARIABLE 

XC=XPNT (LIMP) 

C INITIAL PRFSSURE 

PINL=0. 

C 

SIG = -H. 

IF (XPNT (2) .LT.XPNT(l) ) SIG=-1. 

XWS=XI 

NT=MODE 

NTSV=MODE 

mtp=nt + 1 

NTH-NT/2 

KCMOOO=KCMOOE 

IS = 0 

PI=3. 14169265 

C ************ 

c OIOGNOS T IC FLAGS 

C 

MOIOG = 0 PRINT EVERY DOUBLE OR HALVE OF STEP SIZE 
c MOIOG = 1 NO PRINT 

MD IOG=l 

C 

C MP=0 PRINT EVERY INTEGRATION STEP 




HP=1 NO PRINT 
c MP=1 

C 
C 

c 
c 


kkontm, maximum no integration steps 
kkontm= 300 00 

kprdrv=o print derivatives before profiles 
kprorv=i no print 

<PR0PV=1 

************** 

profile generation points, MUST BE 11 OR 21 

LU-11 

lum=lu-i 

QELZ=0.5/FLOAT<LUM) 

Z w C 1 > = 0 . 3 
00 380 1 = 2, LU 
ZW(I)=ZW(I-1) +OELZ 
380 CONTINUE 

OO 382 1 = 1, NTP 
ZINF<n*0.0 

382 CONTINUE 

00 383 1 = 3 , NTP, 2 
I APG=I-2 

ZINF(I>=. 0C1**IARG 

383 CONTINUE 

00 384 1=1, LU 
ZWSQ=ZW( I) *ZW ( I) 

ZVUV(I, 1)=1. 

ZVW(I,1)=ZW(I) 

00 384 J-3 ,NTP,2 

ZVUVCI, J)=7VUV(I, J-2) * ZWSQ 

ZVWCI , J) = ZVW( I , J-2) *7WSQ 

384 CONTINUE 

00 3 1=1, NTP, 2 
00 3 J= 1 , NTP, 2 

F (I,J)=1.0/(( I+J-l>*2.0**<I*J-l) ) 
as (I, J) = I* (I+J-l) * (I*JM) *2. 0*MH-J+1> 

G(I, J) = (4. 0*1* (I + J + l) +12. 0*( J-l) ) /QS ( I , J) 

M<I,J)*|4.0*IMI«-J + 1> -12.0*1) /QS ( I , J) 

3 CONTINUE 
L C= 0 
KKONT = 0 
IOUTR=0 
0X1=1. E-8*SIG 
2 INFL=0 
ICONT=0 
INFLT=0 
IL IM=0 


KPMAX , PRINT INTERVAL FOR PROFILES ,DI0GN0STIC 
KPNT=0 PRINT PROFILES EVERY KPMAX STEPS 

K PNT = KPMAX -*-1 NO PRINT 
<PMAX = 20 0 
(<PNT = KPMAX+1 



•t q 0 KONT-1 
MC= 0 
MT= 0 

Hi- 0 

M2=0 

M3=0 

M4 = 0 

MCO=0 

KF = 0 

KFF=0 

EF = 1. OE11 

OXM-SIG *0.0025 

DXTE<? = 1 .E-12 

EMAX=1. 0E-04 

EMIN=EMAX/100000. 0 

X= X I 

DX=DXT 

P= PTNL 

08TAIN ORIGINAL VALUES OF DEPENDENT VARIABLES 

DO 4 1=1, NTP 

AC (I) =0 • 0 

8C ( I ) = 0 • 0 

ACD(I)=0. 0 

CONTINUE 

IF (INPRO.NE.2) GO TO 8 
UVPR=UO* VO 

UUUVV=UO*UO*UO+UO*VO*VO 
RI1=59.*UVPR* 0.025/3. 

RI2=59.*UUUW*0.025/3 . 

IF (SIG.GT.O.) PTNL=-RI2/UO 

P=PINL 

IS = i 

XHS=XPNT (2) 

INUM( 1) =1 
PSV(l) =PINL 
RS V ( 1 ) = X I 
ETASV(i> = 0 . 

PTSV(1> =RI2/U0 + PINL 
TQUESV(l) =0. 

8 CONTINUE 
KER = 0 

CALL INPUT ( KER, AUP,R,IFIT , RE , UO , VO , MODE , I NPRO , S T G , p , L U , 

1 KK1,KK2,KK3,PCTU,PCTV) 

IF (KER.E0.2) GO TO 50 
IF (KER.EQ.3) GO TO 51 
00 6 1=1, NTH 
J= I * 2- 1 
K= NTH + I 
YI (I) =AC ( J) 

YI(K)=8C( J) 

6 CONTINUE 
00 16 1=1, NT 
Y ( I > =YI ( I ) 

E ( I ) = 0 « 

15 CONTINUE 



If (INPRO.EQ.2.OR.ILIM.NE.0) GO TO 888 
initial TORQUE ANO DYNAMIC PRESSURE 
C *11*0.0 
* 12=0 • 0 

DO 18 I s 1 ♦ NIP * 2 
DO 18 J=1,NTP,2 

*I1=RI1 + AC (I) *9C( J) /< (I + J-l) *2.Q*MI + J-1) ) 

DO 18 K = 1 » NTP , 2 

rI2=RI2MAC(I> *AC ( J) *AC(K)+AC(I)»BC(J)*BC(K) )/((I+J+K-2> 
1?. 0** (I+J+K-2) ) 

1B CONTINUE 
BB p CONTINUE 
CALL QERV 
CALL ORNT(I) 
kpar=i 

IF (INPRO.EQ.2) KPAR=0 
CALL PROFTL <KPAR,XWS) 

2 0 CALL MOVE 
KKONT=KKONT+l 

IF (KKONT.GT. KKONTM) GO TO 55 
<PNT=KPNT ♦ 1 
IF (MT) 22,21,22 

21 KONT=KONT * 1 

IF ( KONT-4 ) 22,22,23 
21 CALL RKGI 

DO 24 1=1, NT 
E ( I ) =0 . 0 
24 CONTINUE 

IF (MOIOG.EQ.O) CALL PRNTC2) 

GO TO 25 
23 CALL AOMI 

IF (72(1) -EF) 25,25,40 
25 IF (MP.EO.O) CALL PRNT(2> 

IF (ABS(DX) .LT.DXTER) GO TO 32 
C EXAMINE IJ PROFILE FOR INFLECTION 
U I NFL = 0.0 

00 26 I=3 t NTP,2 

UINFL = UINFLMI-1) *AC(I) *ZINF(I) 

26 CONTINUE 

IF (SIG *UINFL) 27,27,28 

27 IF (INFL.EQ.O) GO TO 29 

1 NFL = 0 

ICONT=ICONT+l 
RINF ( ICON T ) = R 
GO TO 29 

28 IF (INFL.EQ.l) GO TO 29 
INFLT=1 
ICONT=ICONT+l 
RINF(ICONT)=R 

I NFL = 1 

29 CONTINUE 

CALL ITRO (MCO) 

IF (MCO. EQ. 1) GO TO 601 

IF ( (ARS(X-XWS) - AB S ( 0 • 99* OX ) ) . GE . 0 • ) GO TO 602 
801 KPAR=1 


GO TO 603 

_ Tf ( K p NT . NE .KPMAX ) GO TO 599 

602 JpNT=0 

, 03 CONTINUE 

6 IF (KPRDRV.EQ.l) GO TO 604 
CALL PRNT ( 2) 

, n4 CALL OROFIL (KPAR,XWS> 
loq CONTINUE 

IF (MCO.NE.O) GO TO 30 
IF (KCMOOE.EQ.O) GO TO 20 
00 =(X-XFIT)*SIG 
IF (OO.LT.O.) GO TO 20 

CHECK ATTEMPT TO REDUCE PROFILE THAT IS INFLECTED 
f IF (CSIG*UINFL) .GT.0. ) GO TO 60 
PEOUCED PROFILE FIT 

L KCN00E=6, 6TH ORDER, KCMODE = 4, 4 TH ORDER 

G KCMM=KCMODE/2-l 

GO TO (700,710), KCMM 
c SWITCH TO 4TH OROER 

700 MOOE-4 

CALL CF IT (IFIT,MODE) 

ILlM=3 
GO TO 720 

c SWITCH TO 6TH OROER 

710 M00E=6 

CALL CFIT ( IFIT ,MODE) 

ILIN=4 

720 CALL PROFIL (0,XWS> 

IS=IS-1 

XFITSV=XFIT 

kcmooe=o 

INPRO-4 

NT=MOOE 

NTP=NT-»-l 

NTH=NT/2 

0X1=1. E-8*STG 

f>INL = P 

XI = X 

00 73 0 1 = 1 , IL IM 
J= IFIT ( I ) 

AUP ( I , 1 ) = ZW( J) 

AUPCI , 2 ) =UW (J) 

AUP( I , 3 ) = VW (J ) 

730 CONTINUE 
GO TO 390 

0 

32 IOUTR=l 
30 CONTINUE 

0 

C PRINT SUMMARY 

C 

COC 6400 TIME RECORO 
ETIME=SECONO( A) 

OFLT=ET IME-BT IME 
CS$$J COC 6400 DATE RECORO 



WHEN=0ATE<0) 

IF (SIG.GT.O.) PRINT 900 
IF (SIG.LT.O.) PRINT 902 
IF (IOUTR.EQ.l) PRINT 925, R 
IF (IUNIFH.EQ.l) INPRO=2 
IF (INPRO.EQ.2) PRINT 919 
IF (INPR0.NE.2) PRINT 920 
GO TO (930,935,940,943), INPRO 
930 PRINT 932, NTSV 
GO TO 945 

935 PRINT 937, NTSV 

IF (KCMODO.EQ.O) GO TO 945 
PRINT 938, KCMOOD , XFITSV 
GO TO 945 

940 PRINT 942, NTSV 
GO TO 945 

943 PRINT 939, NTSV 
PRINT 950 

PRINT 952, <<AUP(I,J), J=i,3), 1=1, LIM) 

PRINT 953 

945 PRINT 921, UO, VO, RE 
PRINT 922 
no 35 1=1, IS 

PRINT 923, PSV(I) ,ET4SV(I) ,°TSV(I) ,PSV(I) ,OPDRSV(I) ,TQUESV(I) , 
1 INUM(I) , VMAG(I) ,VANG(I) 

35 CONTINUE 

IF (INFLT) 36,36,37 

36 PRINT 960 
GO TO 38 

37 PRINT 980 

OO 39 1 = 1 , ICONT , 2 
PRINT 965, RINF(I) 

J»IH 

IF (J.GT. ICONT) GO TO 39 
PRINT 966,PINF(J) 

39 CONTINUE 
PRINT 980 

38 CONTINUE 

PRINT 970, OELT, WHEN 
PRINT 975,KK0NT 

IF (MP.EO.C.OR.MDIOG.EQ.O) CALL PRNT<3> 

GO TO 1 

40 IF(NP.EQ.O) GO TO 452 

IF (HOIOG.NE.O) GO TO 45 
452 IF (LC+NT-58) 450,450,451 

451 CALL PRNT ( 4 ) 

450 PRINT 200 

CALL PRNT ( 2 ) 

PRINT 202 
LC=LC+4 
45 CONTINUE 

IF (K0NT.GT.6) GO TO 55 
0X1=0X172.0 
GO TO 2 
50 PRINT 203 


62 


GO TO 1 
51 PRINT 204 
GO TO 1 

55 MDIOGT*MDIOG 
moiog=q 

CALL PRNT(3) 
mdiog=moiogt 
GO TO 1 
50 PRINT 985 
<CMODE=Q 
GO TO 20 

!01 FORMAT (4I5,F15.0) 

102 FORMAT ( 3F15. 0) 

103 FORMAT (2F15.0) 

104 FORMAT (F15.0) 

135 FORMAT ( 3F15. 0) 

200 FORMAT (35H INITIAL INCREMENT UNSAT ISFACTORY/10X15HEPROR RANGE W 
IAS) 

202 FORMAT < 2X21HHALV I NG THE INCREMENT/) 

203 FORMAT ( 1H1 ,43HSINGUL AR MATRIX ENCOUNTERED, JOB TERMINATED) 

204 FORMAT ( 1 H 1 , 5 X , 54 H MA T CH- UP UNSUCCESSFUL IN 200 F-D STEPS, JOB TERM 
1 INATEO) 

P00 FORMAT < 1H1,5X,13HPUMP ANALYSIS/) 

902 FORMAT ( 1H 1 ,5 X , 16HTUR8I NE ANALYSIS/) 

919 FORMAT (5X , 8 H SUMMARY , /5X , 4 3H INCOMPRESSIBLE FLOW BETWEEN ROTATING 
101 SKS » /5 X , 23H BOYACK INTEGRAL METHOD, 29H MATCHED TO BOYD *-0 START 
1 -UP/) 

920 F ORMAT ( 5 X , 8 H SUMMARY, /5X,43H INCOMPRESSIBLE FLOW BETWEEN ROTATING 
lOISKS,/5X, 23H BOYACK INTEGRAL METHOD,/) 

921 FORMAT <//5X,5H UO -,F8.4,5X,5H VO =,F8.4,5X,5H RE =,F 8 . 4 ,//) 

922 FORMAT (9X ,1HR,6X, 3HET A , 10 X , 2HPT , 12X, 1HP, 10 X , 4HDPDR ,1 1 X , 1MT , 

1 1BX,5HKK0NT,9X, 4HVM AG , 9X , 4HV ANG /) 

92? FORMAT (6X,F5*3,lX,F6.4,lX,F12.6»iX,F12.6,lX,Fi2.6,lX,F12.6,15X,I8, 
1 5X,F8.4,5X,F7.2) 

925 FORMA T(//, 10 X, 34H MINIMUM ALLOWABLE STEP SIZE AT R =,F10.5,//> 

932 FORMAT (5X,I2,17HTH OROER SOLUTI ON ,5 X , 24HPARABOL IC INLET PROFILES/ 
1 ) 

937 FORMAT C5X,I2,17HTH OROER SOL UTI ON ,5X , 22HUNIFORM INLET PROFILES/) 

938 FORMAT ( / 6 X, 2 1H PROFILE REDUCTION TO I1,15HTH ORDER AT R =F6.3/> 

939 FORMAT (5X,I2,i7HTH OROER SOLUTION ,5X , 18HARBITRARY PROFILES/) 

942 FORMAT (5X,I2,17HTH ORDER SOLUTION ,5X , 17HMATSCH ASYMPTOTIC/) 

950 FORMAT ( / , 32X , 1HZ , 10X , 1HU, 10 X , 1H V/ ) 

952 FORMAT C 28 X ,F 8 . 6, F 11 . 6 ,F l 1 . 6 ) 

953 FORMAT (IX,/) 

960 FORMAT (//5X,23HNO U PROFILE INFLECTION/) 

965 FORMAT (6X,24HU PROFILE INFLECTS AT R=F9.6) 

966 FORMAT (6X,26HU PROFILE UNINFLECTS AT RsF9.6> 

970 FORMAT (6X,19HTIME IN EXECUTION =F 8 . 2 , 2X , 7HS ECONOS , 5X , 8 HCOC 6400, 

1 5X,A10/> 

975 FORMAT <6X,19HINTFGRATI0N STEPS =15) 

980 FORMAT <//> 

985 FORMAT ( 1 H 1 ,5 X , 66 HPROFI L E REDUCTION REQUIRED USING INFLECTED PROFI 
1LE, NO REOUCTION ) 

7000 FORMAT (12) 

7005 FORMAT (4F15. 0) 



7015 FORMAT < 3 1 5 , 2 F15. 0 > 
ENO 



SUBROUTINE PROFIL <KPAR,XWS) 

COMMON QS(13,13),F(13,13),G(13,13),H(13 f 13),AC(l3),8C(13>, 

1 Y(12),Y1(12),Y2(12> ,Y3(12> *Y4(12> ,Z(12> ,Z1(12> , 

2 Z2(12) ,Z3(12) ,Z4(12),YI(12) ,E(12) ,XI,XC,X,X1,X2, 

3 X 3 » X4, DXI,DX,DXM,EMIN, EMAX , KONT , LC , MP , MC , Ml , M2 , 

4 M3,M4,MT,NT,RE,U0,V0,DPDR,P,P1,P2,P3,P4,PG1,PG2, 

5 MOOE,NOPNT,INPRO,KKONT,MDIOG,KF,KFF,KPNT,SIG, ALPRM, 

6 ATES(64,65) ,XTES(64) , 

7 RSV (31) ,ETASV(31),PTSV(31) ,PSV(31) , TQUESV(31) , 

8 UPOPSV (31) , ACD(13) ,ZINF(13) , ZVUV (21 , 13) , ZVW (21 , 13) , 

9 I NUM (31 ) »RINF(11) , AUP ( 7 , 3 ) , UW ( 21) , VW ( 21 ) , RI1 , 

1 ZW(21) ,WW(21> ,Q(13,13> ,XPNT(31) , RI2 , I S , L U , NTP , P I , NTH , 

2 VMAG(31> ,VANG( 31) 

EQUIVALENCE <X,R> 

RI 3-0 . 0 

IF (SIG .GE.O.) R 1 2= 0 • 0 
00 610 I-1,NTP,2 
DO 610 J-1,NTP,2 
RI 3=RI 3+ AC (I)*BC(J)*F(I,J) 

IF (SIG • L T • 0 • ) GO TO 610 
DO 605 K=1,NTP,2 
IJK=T+J+K-2 

RI 2=R I2+ ( AC (I ) *AC(J> *AC(K) +AC(I) *8C ( J) *BC ( K) ) / ( ( I JK) * 

1 2 • 0**1 JK) 

60S CONTINUE 
610 CONTINUE 

IORQ=-4.0*PI*(RI1-RI3*R*R) 

PTOTL=RI2/UO 

IF (SIG .GE.O.) PTOTL=PTOTL*R 
PTOTL=PTOTL+P 

ETA=-TORQ/ (2.*PI*UO*PTOTL) 

IF (SIG .GE.O.) ETA=-1./ETA 
C STORE PERFORMANCE PARAMETERS 

IF (KPAR.NE.l) GO TO 620 
KP AR= 3 
IS=IS+1 

C NEXT X PRINT STATION 

ISS=IS+1 
XWS=XPNT (ISS) 

INUM( IS) =KKONT 
PSV(IS) 

OPD p SV (IS) =DPDR 

RSV(IS) =R 

ET ASV ( I S) =ETA 

PTSV( IS) =PTOTL 

TQUESV(IS)=TORQ 

CALL VECTOR(VELMAG,VELANG) 

VMAG(IS)=VELMAG 
V A NG ( IS) = V ELA NG 
620 CONTINUE 

IF (NOPNT.EQ.l) RETURN 
IF (SIG.GT.O.) PRINT 900 
IF (SIG.LT.O.) PRINT 902 
PRINT 901 

GO TO (830,835,840,845), INPRO 



g30 PRINT 932, NT 
GO TO 850 
0 35 PRINT 937, NT 
GO TO 850 
040 PRINT 942, NT 
GO TO 850 
845 PRINT 939, NT 
050 CONTINUE 

PRINT 909, UO , R,VO,OX , RE,KKONT 
PRINT 910, OPDR,P,PTOTL,ETA,TORQ 
PRINT 904 
SUM- 0 . 

OO 505 1 = 3, NT P, 2 
IARG=I-3 

505 SUM=SUN+ (FLOAT ( I- 1 ) *FL OA T ( 1-2 ) * AC ( I ) ) /<2**IARG) 

SUM=SUM/RF+R 

C0MPT=ABS< (SUH-OPOR) /DPOR) *100. 

OO 506 J = 1 , NT H 
K=2*J-1 
ACO(K)=Z <J) 

506 CONTINUE 

IF (NT.EQ.4) ACD<5)=-16.MZ(l)+Z(2)/4.) 

IF (NT.EQ.6) ACD(7> =-64. * < Z(1 > +Z (2) /4.*Z< 3) /16. ) 

IF (NT.EQ.6) AC0(9) = -256.*(Z(l)>Z(2)/4.4>Z(3)/16. + Z(4)/64.) 

IF (NT.GT. 8) ACD(NTP)=ALPRH 
OO 510 1=1, NT P, 2 
OO 510 J= 1 , NTP , 2 
FJM=FLOAT ( J-l) /FLOAT (I) 

Q(I,J)=AC(I)* ACO( J> -FJM*AC (J) MAC (I) /R*ACO<I) ) -BC(I) *BC(J)/R 
510 CONTINUE 
CSUM=Q • 0 

OO 515 1=1, NTP, 2 
OO 515 J = 1 , NTP , 2 

CSUH=CSUM+FLOAT(I) *Q(I,J)/QS(I, J> 

515 CONTINUE 

CSUM=-12. *UO/ (R*RE> -24.*CSUM 
COMPTC= APS ( (C SUM- DPOR) /OPOR) *100. 

UW (1) =AC (1 ) 

VW ( 1) =9C ( 1 ) 

WW(1)=0.0 

OO 521 1 = 2, LU 

UH(I)=0. 0 

VW(T) =0 . 0 

WW(I)=0. 3 

OO 518 J= 1 ,NTP , 2 

UW(I)=UH(I)+AC(J) *ZVUV(I, J) 

VW<I) = VW ( I ) +9C ( J) *ZVUV(I,J> 

FJ=J 

WW(I)=WH(I)-(AC(J) /RF ACO(J) ) * ZVW ( I , J) /FJ 
518 CONTINUE 
521 CONTINUE 

OO 525 1=1, LU 

PRINT 905, ZW<I) ,UW(I) ,VH(I) ,WW(I) 

525 CONTINUE 
PRINT 9Q6 



00 530 1=1* NT P, 2 
K=I-1 

PRINT 907* K* AC ( I ) ,BC CD 
530 CONTINUE 
PRINT 915 

PRINT 911* COHPT * CONPTC 

IF (MP.EQ.O.OR.HOIOG.EQ.O) CALL PRNT ( 4 ) 

RETURN 

900 FORMAT < 1H 1 ,5 X , 13HPUMP ANALYSIS/) 

901 FORMAT < 1 8X , 32H REGULAR INTEGRAL SOLUTION , / / , 14X , 43HFULL ADMI 
1SSION, INCOMPRESSIBLE LAMINAR FLOW) 

902 FORMAT ( 1H1 ,5X , 16HTURBINE ANALYSIS/) 

904 FORMAT ( 13 X,1HZ, 11 X , 1HU, 13X , 1HV, 13X ,1HW*/) 

905 FORMAT ( 10 X ,F6 . 3 , 4 X , F 1 0 . 6 , 4 X , F 1 0 . 6 , 4X , F 10 . 6 ) 

906 FORMAT ( /// »15X ,35HU ANO V POLYNOMIAL COEFFICIENTS ARE , / / * 16X , 1H I , 
19 X * 4HA ( I > ,12X,4HB(I) ,/, 15 X , 3H ***, 7X , 6H**** **, 1 OX , 6H* *♦***,// ) 

907 F0RMAT(15X,I2,4X,E12.5,5X,E12.5> 

909 F0RMAT(///,16X,5HU0 = ,F 7 . 3 , 1 IX , 6HR = , El 2 • 5 , /16X , 5H VO =,F7«3, 

1 1 1 X » 6H0R =,E12.5*/16X,5HRE = , F7 . 3 , 11X, 6HSTEP =,H2) 

910 F0RMAT(//*26X* 6 HOP OR = , F 10 . 5 , /26X , 6HP = , F10 . 5 } /26X , 6HPT = , 

1F10.5,//26X,6HETA -F 10 . 5 , /26 X , 6HT =,F10.5,/> 

911 FORMAT C / , 15 X » 12HC/L -TO- WALL=F8 . 2 , 2 X , 7HPC OIFF , / , 15X , 1 2HC /L- TO-Fl 
10W=F8.2,2X,7HPC OIFF) 

915 FORMAT ( / , 15X , 31H PRESSURE DERIVATIVE COMPARISONS/) 

932 FORMAT (5X,I2,17HTH ORDER SOLUTION, 5X ,24HPARAB0LIC INLET PROFILES/ 
1) 

937 FORMAT (5X,I2,17HTH ORDER SOL UTION , 5 X , 22HUNIFORM INLET PROFILES/) 
942 FORMAT (5X,I2,17HTH ORDER SOLUTION, 5X ,17HMATSCH ASYMPTOTIC/) 

939 FORMAT (5X,I2,17HTH ORDER SOLUTION ,5X , 1SHAR8ITRARY PROFILES/) 

END 



SUBROUTINE INPUT ( KER , AUP , R , IF I T , RE , UO , VO , MOOE , INPRO , S IG , P , L U , 
1 KK1,KK2,KK3,PCTU,PCTV) 

£ THIS ROUTINE COMPUTES THE INPUT U,V VELOCITY PROFILES 

C COMMON FILL(889),A(64,65> ,X<64) ,0UM(212) ,ZVUV(21,13> 

OIMENSION AUP(7,3) ,AC(13) ,BC(13) ,ZT(2l),UT(2i) ,VT(21) 

OIMENSI ON 72(7) ,74(7) ,Z6<7) ,Z8(7) ,710(7),Z12(7) ,Z14(7) 
OIMENSION B(7,7),IFIT(7) ,VC(21) 

EQUIVALENCE ( FILL (677) , AC ( 1) > , (FILL (690) ,BC ( 1 ) ) 

EQUIVALENCE ( OUM( 156) , PORI V) 

NTP=MODE+l 

LlM=MODE/2+l 

LIMP1=LIM+1 

GO TO (1,2,3, 4) , INPRO 

C 

c PARABOLIC VELOCITY PROFILES 

1 UMAX=1.5*U0 
VMAX=1.5*V0 
AC ( 1) = UM AX 

AC (3)=-4. *UMAX 
BC (1> =VMAX 
BC (3)=-4.*VMAX 
300 RETURN 
C 

r UNIFORM PROFILES 

2 CONTINUE 
KER=0 
K1=KK1 
K2-KK2 

50 CALL TESLA (UO , VO , RE , S IG , KER, ZT , U T , VT , P , UA , R , K1 , K2 , KK 3 , PORI V ) 
IF (KER.GE.2) RETURN 
UTEMP-UO 
DO 100 1=1, LIM 
J=IFIT (I) 

AUP (I, 1)— ZT (J) 

AUP(I ,2) =UT(J) 

AUP(I,3)=VT(J) 

19C CONTINUE 

on 113 1=1, LTM 
A ( I , L I MPl ) =AUP (1 , 2) 

7 . 7 = AUP (1,1) 

11 SQ=7Z*7Z 
0(1 , 1 ) = 1 • 

A(I,1)=1« 

DO 110 J=?,LIM 
JM 1= J-l 

1 ( T , J ) = B ( I , JM 1 ) *Z ZSQ 
A(I,J)=B(I ,J) 

110 CONTINUE 
GO TO 400 
C 

C MATSCH 3RD APPROX ASYMPTOTIC PROFILES 

3 0ELTAZ=0.5/FL0AT (LIM-1) 

DO 210 J=1,LIM 



F J- J 

ZZ =0ELTAZ*(FJ-1.) 

RE2=RE*RE 

RE3=RE*RE*RE 

UR=U07R 

U2R=tU0*U0)/(R*R*R) 

U3R= <UO*UO*UO) 7<R*R*R*R*R> 

U4R=(UO*UO*UO*UO> 7<R*R*R*R*R*R*R) 

Z2(J)=ZZ*Z7 

Z4(J)=Z2(J)*Z2<J> 

Z6(J)=Z4< J)*Z2(J> 

Z8(J)=Z4(J)*Z4(J) 

Z10(J)=Z6( J)*Z4(J) 

Z12( J) = Z6 ( J) *Z6 < J) 

Z14(J> =Z8 ( J)*Z6(J> 

T1 = -12.MZ2(J> *.5-.125>*UR 

T2=-RE*U2R* (1.2*Z6(J> -1. 5*Z4 ( J) f 99 . *72 l J) /2 8 0 . -3 . /224 . ) 

T3 = -RE2*URM-Z6<J> 715 •♦•25*Z4 ( J) -39.*Z2(J> /560 . + 19. /6 7 20 . > 

T4 = RE3*U2RM-Z10< J> 790 . *3 . *Z8 ( J) 756.-23. *Z6 (J) /240 . *5 . * Z4 ( J ) 7 64 . 
i -29793. *Z2< J)/17740 80. +4361. /7 096320.) 

T5=RE2*U3RM-48.*Z10< J) 7175. +9 . *Z 8 ( J) / 14. -171. *Z6<J> 7 3 50. 
i ♦99.*Z4(J)7560.-2217.*Z2(J) 786 240 . +5301 . 768992 0 0 . > 
T6=RE3*U4RM-162.*Z14(J> 715925. + 129.* Z12(J>7385 0.-409. *Z10(J>7 
1 9800. + 893.* Z8(J> 739200.- 879. *Z6 (J) 7156800 . +99. *Z4< J ) 7 12 5440 . > 

T7=RE3*U4RM-876231.*Z2(J) / 1 2 5565540 0 0 . + 96 151 . /5 0 2261 76 00 3 . ) 

UA SYN = T 1+T2FT3+T4+T5+T6+T7 

VASYH= R-(RE*UR*( Z4(J) -1.5*Z2 (J) ♦ (5./16.M ) -<RE2*U2R* ( (3. *Z8 ( J) / 

1 70.)-.l*Z6(J)+(33.*Z4(J)/560.)-(3.*Z2(J) 7224. )+< 19. 717921. > ) ) 

2 ♦ <RE3*U3RM(4.*Z12(J)7385.)-(4.*Z10< J) 7105.) +(12. *Z8( J) 7245.)- 

3 ( 57. *Z6 < J) /2100. > +<3.*Z4< J> / 448. >-(3127.720697600.) ) ) 

A ( J,LTMP1> = UA S YM 

AUP(J,3) =VASYM 

no 210 1 = 1 , LIH 

K= 1*2-2 

B< J,I>=ZZ**K 

A(J,I)=B(J,I) 

210 CONTINUE 
GO TO 400 
C 

C AR8ITRARY VELOCITY PROFILES OR AC ( I ) ANO BC(I> ON AUP(I,J> OATA CARDS 

4 IF « AUP ( 1 , 1 ) . G T , 0 . 5) GO TO 500 
00 310 1 = 1. LIM 
ZZ=AU°< I. 1) 

A { I ,LIMP1) =AUP (I, 2) 

Z7SQ=ZZ*ZZ 

9(I,1)=1. 

A ( 1,1) =1 , 

00 310 J= 2 , LI M 
JH 1= J- 1 

B(I,J)=8(I,JM1)»ZZSQ 
A (I, J) =B (I , J) 

310 CONTINUE 

C 

400 CALL GAUSSN < A , X , 64, 65 , LI M , KER ) 

IF (KER.EQ.2) RETURN 



00 410 1=1, LIM 
K=I*2-1 
AC <K) = X ( I) 

410 CONTINUE 

DO 420 1 = 1, LIM 
A(I,LIMP1)=AUP(I,3> 

00 420 J = 1 , LI M 
A (I,J)=B(I ,J> 

420 CONTINUE 

CALL GAUSSN ( A , X , 64, 65 , L IH, KER) 

IF (KER . EQ . 2) RETURN 
DO 430 1 = 1 , LIM 
K= 1*2-1 
90 (K> =X ( I ) 

430 CONTINUE 

C COMPUTE UO BY INTEGRATION OF U PROFILE 

UO = 0. 

DO 44 G 1 = 1, LIM 
K= 1*2-1 
FK=2**K 
FKK=FK*K 
UO=UO+AC (K) /FKK 
440 CONTINUE 
UO=UO*2.*R 

I F ( INPRO • E0.2 ) GO TO 600 
RETURN 

C 

500 UO=0. 

00 510 1=1, LIM 
K=2*I-1 

AC ( K) = AUP ( I ,2 ) 

8C(K)=AUP(I,3) 

FK=2**K 
F K K= F K * K 
UO=UO+AC(K) /FKK 
510 CONTINUE 
UO = UO * 2 • *R 
RETURN 
600 CONTINUE 

00 610 1 = 1 , LU 
VC (I) = 0. 

00 610 J = 1 , NT P , 2 

VC (I) =VC (I ) +BC ( J) * ZVUV < I , J > 

610 CONTINUE 
M A R= 1 

IF (LU. EQ, 11) MAR=2 
J=0 
K V= 0 

00 620 1 = 1, 19, MAR 
J=J + 1 

VOUM= ABS ((VC(J)-VT(I))/VT(I))*100. 

IF (VOUM.GT.PCTV) KV=KV+1 
620 CONTINUE 
KU= 0 

UOIF=ABS ( (UO-UA) /UO >*100. 



/ u 


IF (UOIF.GT.PCTU) KU=1 

IF (KK3.EQ.1) PRINT 630, ( I , AC C I ) , I, BC (I) , I=1,NTP,2> 

6 30 FORMAT ( 5X , 3H AC ( , 12, 2HI *E18 • 1 1 , 5X , 3HBC < » 12 , 2H> =E1«. H > 
IF (KU.EQ. O.AND.KV.EQ. 0) RETURN 
K1=K1+10 
KER=3 

IF ( K 1 • GT .20 0 ) RETURN 
<2 = 10 
KER=-1 
UO=UTEMP 
GO TO 50 
END 



SUBROUTINE TESLA ( UO , VO , RE , SIG, KER, Z , U, V, P , U A ,R , KK1 , KK 2 , KK3 , PORI V ) 
COMMON FILL(889),A(64,65> ,X<64> 

OIHENSION U<21), V<21>, W(21>, OU(2i), OV(21), OW<21> 

(DIMENSION S (21 ) ,Z(21) 

KCMAX=KK1 

IF <KK2 .LE.O. 0R.KK2.GT.KK1) KK2-KK1 

KCCMX=KK2 

KPRNT=KK3 

IF (KER.EO.-l) GO TO 350 
FJW=21 • 

JW=2l 

OZ = .5/ (F JW- 1 . ) 

JWM4 - JW-4 
JWM3 - JW-3 
JWM2 = JW-2 
JWM1 = JW-1 
JWP1 - JW+1 
JWP2 = JW*2 
JWP3 = JW-t-3 
JW2M6 = 2* JW- 6 
JW2H5 = 2* JW- 5 
JW2M4 = 2 * J W- 4 
JW2N3 = 2* JW-3 
JW2 = 2* JW 
JW2P1 - 2*JW+1 
JW2P2 = 2*JW+2 
JW2P3 = 2*JW+3 
JW3N6 = 3*JW-6 
JW3M5 = 3 * JW- 5 
JW3M4 = 3* JW-4 
JW3M3 = 3* JW-3 
JW’*2 = 3* JW- 2 
JW3M1 = 3* JW- 1 
JW3 = 3* JW 
JW3P1 = 3* JW-*- 1 
JW3P2 = 3* JW+ 2 
JWM1D2= JWM 1/2 
R = 1. 

C UNIFORM VELOCITY PROFILES AT R=1 
00 65 J=1,JW 
U(J) - UO 
V(J) = VO 
W(J> - o. 

65 CONTINUE 
U(JWJ = 0. 

VCJW) = R 

l INTEGRAL (USQ*-VSQ) at r=i, z=o TO Z= • 5 

51 00 52 J=l, JW 

SCJ) - (U(J)*U(J>) + (V(J)*V(J)> 

52 CONTINUE 
x° = 0. 

00 56 J=l, JWM102 
JJ = 2* J 
JJJ = JJ*-1 

XPT = (4.*S(JJ)» ♦ ( 2 • *S ( JJJ ) ) 



XP = XP+XPT 
56 CONTINUE 

XP = (XP+S ( 1) -S CJM) ) * DZ / 3. 

TERM1 = 2. * XP 

TERM2 = -2. * (U(JWM2) - 4. * U(JWMl) + 3. * U(JW)) 
OPOR = TERM1 + TERM? 

PDRIV=OPOR 
2 0 C CONTINUE 
ETA = 0. 

T - 0. 

205 DO 210 J=1,JW 
FJ = J 

Z< J) = <FJ-1. ) * OZ 
210 CONTINUE 

C TORQUE, UV PROOUCT, R=1 

C INTEGRAL UV AT R= 1 , Z=0 TO Z=.5 

250 XO = 0. 

00 255 J=i, JHM102 
JJ = 2*J 

JJJ = JJ + 1 

XOT = (4 . *U(JJ) *V ( JJ) +<2.*U(JJJ> *V(JJJ)) ) 

XO = XO + XOT 
255 CONTINUE 

XO = (XO+ ( U(1 ) *V( 1) ) - <U< JW) *V ( JH) ) ) ♦ OZ/3 . 

C DYNAMIC PRESSURE, R=1 

C INTEGRAL (USQ+VSQ)*U AT R=l, Z=0 TO Z=.5 
DO 260 J = l, JH 

SCJ) = ( (U(J> *U(J) ) + < V(J) *V< J) ) ) *U(J) 

260 CONTINUE 
XP = 0. 

DO 265 J=1 , JWM102 
JJ = 2* J 
JJJ = JJ + 1 

XPT = <4.*S(JJ))+(2.*S<JJJ)) 

XP = XP+XPT 
265 CONTINUE 

XP = ( XP + S ( 1) - S( JH) ) *OZ/3. 

PT = P + (XP/UO) 

01 = RE*OZ *DZ 
8D1=1./81 
002=2. /B1 

82 = ?.*OZ 
020=1. /R2 
0R=1.E-6*SIG 
KCC = 0 
KC = 0 
GO TO 1000 
350 R = R + DP 

ROR=l./OR+l./R 
OPORP = DPOR 
DO 365 M= 1 , JH3P1 
00 360 N = 1 , JW3P2 
A ( M,N) = 0. 

360 CONTINUE 
365 CONTINUE 



00 370 M - 1, JWM2 
A ( Mf M+ 1 ) = <U <M+1>/DR) ♦ BD2 
A(M»M> * -(W(Mfl)/82) - 801 

A(M v M + 2> = <W<M*1>/B2> - 001 
MPJWP1=M+JWP1 

A(M,MPJWP1) = -2. * V(H+1» / R 
MJW2P1=M+JW2P1 

A(W,HJW2P1 ) = (U(H4-2)/92) - <U(M)/B2> 

A ( M » JW3P 1) = 1. 

A (M,JW3 P2> = < V(M + 1> *V(M + 1>/R) - (W(H+i> *U(M+2) /82) -MW< M + i) *U(M) /R2) 
1 ♦ ( (U(M+2) +U(M) > /Bl) - ( <2.*U<M*i> ) /Bl) 

370 CONTINUE 

00 380 M = JW Ml y JW2M4 
MM JWM3=M- JWM3 
A ( M , HMJ WM3 ) = V (MMJWM3) /R 

A ( M » M+- 3 > = <U(MMJWM3> /DR) ♦ ( U(MMJWM3 ) /R) ♦ 802 

A(M,M+2) = -( W (MMJWM3) /B2) - B01 

A(M,M+4) = (W (MMJWM3) /B2) - B01 

MMJWM4=M-JHM4 

MMJMM2=M-JWM2 

NPJWP3=M+ JWP3 

A<M,MPJW®3) = ( V ( MMJWM4) - V ( MM JWM2) ) /B2 

A ( M» JH3P2) = ( - (U ( MMJ WM3 ) * V ( MM JWM3 ) ) /R) - ( W ( MMJHM3 ) MV (MMJWM4) 

1 - V ( MM JWM2) )/B2)+ < V (MMJWM4) /Bl) - ( 2**V (MMJWM3) /Bl) ♦ ( V ( MM J WM2 > / 81 ) 
380 CONTINUE 

00 390 M=JW2M3, JW3M6 
MJW2M5=M- JW2M5 
A ( M, MJW2M5 ) = ROR 
A ( M,M+6) =820 
A(M,M+4) =-820 
M JW2M6=M- JW2M6 
MJW2M4=M-JW2M4 

A(M,JW3P2) = (-UCMJW2M5 > /P) - ( < W { M JW2M6 )-W(MJW2M4 ))/82) 

390 CONTINUE 

A ( JW3M5 » JW ) = 1. 

A C JW3M4 ♦ J M2) = 1. 

A ( JW3M4, JW3P2) = R-V(JW) 

A ( JW3M3 t JM 3 ) = 1. 

A(JW3M2,1) = -3. 

A(JW3M2,2> = 4. 

A ( JW3M2 » 3) = -1. 

A ( JW3M2, JW3P2) = 3**U(1)-4»*U(2)4-U(3) 

A( JW3M1, JWP1) = -3 
A( JW3M1, JWP2) s 4. 

A ( JW3M1, JWP3) = -1. 

A( JH3M1 , JW3P2) = 3.*V(1)-4.*V(2)+V(3) 

A ( JW3 » J W 2P 1 ) = 1. 

A ( JW3 * J W3P2) = -W(l> 

A ( JW3P1, JW3M2) = 1. 

A( JW3P1, JW3M1) = -4. 

A(JW3P1, JW3P2) = 4.*W ( JWM1) -W ( JWM2) 

CALL GAUSSN <A t X, JW3P1, JW3P2, JH3P1, KER) 

IF (KER* EQ .2) RETURN 
395 00 400 J=1,JM 
0U(J) = X(J) 



JPJW=J+ JW 
JPJW2=J+JW2 
OV(J) = X(JPJW) 

OH(J) = X(JPJH2) 

OPOR = X (JW3P1) 

400 CONTINUE 

00 410 J=1,JW 

U(J) = U(J) ♦ OUCJ) 

V(J) = V(J) + QV(J) 

W(J) = W(J) + 0W(J) 

410 CONTINUE 
U(JW) = 0. 

V(JW) = R 
W(JW) = 0. 

CALCULATE PRESSURE 

500 P = p ♦ (DPORP+OPOR) * .5 * DR*SIG 
KCC=KCC+1 
<C=KC+1 

IF (KC.GT.101) GO TO 520 
FKC=KC 

0R=(l.E-6+9.99E-8*FKC*FKC> *SIG 
GO TO 530 
520 DR=.001*SIG 
530 CONTINUE 

IF (KC.GE. KCMAX) GO TO 600 

599 IF (KCC.NE.KCCMX) GO TO 350 
KCC=0 

600 XI = 0. 

IF (SIG.GT.0.) XP=0. 

00 610 J=l, JWM102 
JJ = 2»J 
JJJ = JJM 
TORQUE, R 

XI T = ( 4 • * U ( J J ) *V < JJ) ) ♦ (2.*U( JJJ) *V ( JJJ) ) 

XI = XI+XIT 

IF (SIG.LT.O.) GO TO 610 
DYNAMIC PRESSURE, R 

XPT=4.» (U( JJ) *U(JJ)*U(JJ) +V( JJ) *V (JJ) »U(JJ) )+2 .* (U(JJJ) *U< JJJ) 
1 U(JJJ)+V(JJJ)*V<JJJ)*U(JJJ)> 

XP=XP+XPT 

610 CONTINUE 

XI = (XI+ (U(l) *V(i) ) - (U(JH) ( JW) ) ) * DZ/3. 

IF (SIG.LT.O.) GO TO 611 

XP=(XP+ U(1)*U(1)*U(1)+V(1)*V(1)*U(1)-U(JW)*U(JW)*U(JW>- 
1 V(JW)*V(JW)*U(JW> )*0Z/3. 

611 CONTINUE 

PT = P*-XP*R/UO 

IF (SIG.LT.O.) PT=P+XP/UO 
RSQ = R*R 

T = -(XO- (PSQ*XI) ) * 12.566370606 
WI - 6.283185307 * PT * UO 
ETA=WI/T 

IF (SIG.GT.0. .AND. R.EQ.l.) ETA=0. 

IF (SIG.LT.O.) ETA=-T/WI 
IF (KPRNT.EQ.l) GO TO 1001 



/ b 


GO TO 1 D 8 5 

iOOO IT" (KPRNT.EQ.O) GO TO 350 
tOOl °RINT 1 

1 FORMAT (1H1) 

PRINT 1005 

1 0 05 FORMAT (50H FLOW BETWEEN ROTATING DISKS, INCOMPRESSIBLE FLUID,/) 
PRINT 1010 

tC 10 FORMAT ( 34H BOYD ANALYSIS, PROGRAMMED BY RICE,/) 

1020 PRINT 1021 

1021 FORMAT (26H UNIFORM VELOCITY PROFILES,/) 

1040 PRINT 1041, UO, VO, RE 

1041 FORMAT (IX, 4HU0 = F8 . 4 , 5 X , 4H VO =F8.4,5X, 4HRE =F8.4,/) 

PRINT 1043, R , DR ,KC 

1043 FORMAT (4H R =F12 . 8 , 1 0 X , 4H0R =F 1 2 . 8 , 9 X , 3HKC -1 5 / ) 

PRINT 1045, DPOR, P, PT 

1045 FORMAT ( 7H OPDR = F 15 . 6 , 5 X , 3HP =F15. 8 , 5X,4HPT =F15.8,/) 

PRINT 1050, T, ETA 

1050 FORMAT (4H T =F15 . 8 , 1 OX , 5HET A =F10.5,//) 

PRINT 1055 

1055 FORMAT (3X,2H Z,14X,2H U,17X,2H V,17X,1HW,/) 

DO 1070 J- 1, J W 

PRINT 1060, Z(J), U(J), V(J), WCJ) 

1060 FORMAT ( F6 . 3, 3 ( 4X , F15 . 8 ) ) 

1070 CONTINUE 
1085 CONTINUE 
C COMPUTE CONTINUITY 
UA = 0. 

DO 1300 J=l,10 
J J = J* 2 
JJJ=JJ+1 

UAA=4. *U( JJ)+2.*U (JJJ) 

UA = UA + UA A 
1300 CONTINUE 

UA-(UA+U (1) -U (21) ) *2. *R*Q. 025/3. 

COMPR=( (UA-UO) /UO) *100. 

IF (KPRNT.EQ.0) GO TO 1090 
PRINT 1080 ,UA,COMPR 

1080 FORMAT (6X,15HF-D CONT INUIT Y=F 1 0 . 6 , 5 X , 5H0IFF-F6 . 3 , 3H PC) 

1090 IF (KC.NE. KCMAX) GO TO 350 
RETURN 
ENO 



SUBROUTINE CFIT ( IFIT ,MOOE) 

CUPVE FIT POINTS CORPESPONTING TO OROER OF FIT 
OIHENSI ON TFI T (7) 

NOTE THIS SUB MUST HAVE 11 OR 21 PT FIT ADDED 

<0IR=M0DE/2-l 

GO T0(4, 6, 8,10,12) , KDIR 

4TH ORDER INTEGRAL REDUCTION (11 PT) 

IF I T ( 1 ) =1 
IFIT(2>=7 
IFIT( 31 =11 
RETURN 

6TH OROER INTEGRAL REDUCTION (11 PT) 

IF IT ( 1 ) = 1 
IFIT (2) =5 
IFIT ( 3 ) = S 
IFIT<4) =11 
RETURN 

8TH ORDER F-D MATCH-UP (21 PT) 

IFIT(l) =1 

IFIT(2) =6 

IFIT(3)=11 

IFIT ( 4) =16 

IFIT (5)=21 

RETURN 

CONTINUE 

RETURN 

CONTINUE 

RETURN 

END 
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SUBROUTINE GAUSSN ( A , X , NN , MM , N , KER ) 

DIMENSION A (NN,MM) ,X(NN) 

A IS AN N BY N + l MATRIX OF COEFFICIENTS FOR A SET OF N 
SIMULTANEOUS LINEAR ALGEBRAIC EQUATIONS. THE CONSTANTS ARE 
STORED IN COLUMN N+l. X IS THE N DIMENSIONAL ANSWER VECTOR 
THIS SUBROUTINE HAS VARIABLE DIMENSIONS (NN AND MM) 

A IS DIMENSIONED A(NN,MM) AND X IS DIMENSIONED X(NN) 

NPM-N+1 

EP-l.F-27 

10 DO 34 L = 1,N 
KP=0 

Z=0.0 

DO 12 K = L , N 

IF(Z-ABS( A (K ,L) ) ) 11, 12, 12 

11 Z= ABS ( A ( K , L ) ) 

KP=K 

12 CONTINUE 
IF(L-KP) 13,20,20 

13 00 14 J-L,NPM 
Z=A(L, J> 

A<L,J)=A(KP, J) 

14 A(KP,J)=Z 

20 IF(ABS(A(L,L) )-EP)50,50,30 

30 IF(L-N) 31,40,40 

31 LP1=L+1 

DO 34 K = L PI , N 
IF(A(K,L) >32,34,32 

32 RATI0=A(K,L) / A (L , L ) 

DO 33 J=LP1,NPM 

33 A <K, J> =A (K,J> -RATIO+A (L , J> 

34 CONTINUE 

40 DO 43 1 = 1, N 
II=N+1-I 
S=0. 0 

IF(II-N) 41,43,43 

41 IIP1=II+1 

OO 42 K= I IP1 , N 

42 S=S+A ( II, K)*X (K) 

43 X (II> = < A(II,NPM)-S)/A (II, II) 

KER = 1 
GO TO 75 
50 KER- 2 
75 CONTINUE 
RETURN 
END 
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SUBROUTINE ITRO (<) 

THIS ROUTINE CHECKS THE CURRENT VALUE OF X WITH THE SPECIFIED 
CUT-OFF VALUE OF X 

COMMON QS (13, 13) ,F (13,13) ,G( 13,13) ,H( 13 , 13 ) , AC (13) ,BC (13) , 

1 Y ( 12), Y 1(1 2) ,T2(12) ,Y3(12> , Y4 ( 12) , Z (12) ,Z1(12) , 

2 Z2(12) ,Z3( 12) ,Z4(12> , YK12) ,E(12) , X I , XC , X , XI , X 2 , 

3 X3, X4,0XI,DX,DXM,EMIN,EMAX,K0NT,LC,MP,MC,M1,M2, 

4 M3,M4,MT,NT,RE,U0,V0,DPDR,P,P1,P2,P3,P4,PG1,PG2, 

5 MOOE,NOPNT, INPRO,KKONT ,MOIOG,KF,KFF,KPNT,SIG, ALPRM 
EQUIVALENCE (X,R) 

IF (SIG MX+DX-XC)) 2,1,1 

1 IF (M4) 4,5,4 

2 K= 0 

3 RETURN 

5 DX-XC-X 
KONT =1 
M4=l 
GO TO 3 

4 K= 1 

GO TO 3 
END 
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SUBROUTINE move 


THIS ROUTINE RE-STORES THE VARIABLES OF PREVIOUS STEPS BEFORE 
PROCEOING TO CALCULATIONS FOR THE NEXT STEP 

COMMON QS <13, 13) ,F (13,13) ,G< 13,13) ,H< 13,13 ), ACC 13) ,BC (13) , 
Y(12), Yl(12> ,Y2(12) ,Y3(12) ,Y4(12) ,Z(12) ,Z1<12> , 

Z2(12) ,Z3(12) ,Z4(12) ,YI(12) ,E(12) , XI , XC , X , XI, X 2 , 
X3,X4,0XI,0X,DXM,EMIN,EMAX,K0NT,LC,MP,MC,M1,H2, 
M3,M4,MT,NT,RE,U0,V0,DPDR,P,P1,P2,P3,P4,PG1,PG2, 

mode,nopnt,inpro,kkont,mdiog,kf,kff,kpnt,sig, alprm 

EQUIVALENCE ( X , R) 

P4=t»3 
P3-P2 
P2 = °l 
P1 = P 
PG2=PG1 
PGl^OROR 
X4 = X3 
X3 = X2 
X2 = X1 
X 1 = x 

DO 1 1=1, NT 
Y4 ( I ) = Y 3 ( I ) 

Y3 (I) =Y2 (I) 

Y2(I)=Y1(I) 

Yl(II-YCI) 

Z4(I)=Z3(I) 

Z3(I)=Z2 (I) 

Z2(I)=Z1(I) 

Z1(I)=Z(I) 

1 CONTINUE 
RETURN 
ENO 



SUBROUTINE PRNT <K> 


' THIS ROUTINE PROVIOES PRINTEO OUTPUT FOR EACH INTEGRATION STEP 

' WHEN THE INPUT VARIABLE MP=0. THIS PRINTOUT MAY BE CONSIDERED 

l r LARGELY DIAGNOSTIC IN NATURE AS IT DOES NOT INCLUOE U, V, W OR 
r TURBINE PERFORMANCE INFORMATION 

l j 

C COMMON US<13, 13) ,F (13, 13) ,G(13,13) ,H(13,13) ,AC(13> ,BC (13) , 

1 Y(12),Y1(12),Y2<12>,Y3<12>,Y4(12),Z(12),Z1<12), 

2 Z2(12) ,Z3(12) , Z4<12> ,YI(12) ,E(12) , XI , XC, X, XI, X2, 

3 X3,X4,DXI,DX,DXM,EMIN,EMAX,K0NT,LC,MP,MC,M1,M2, 

4 M3,M4,MT ,NT ,RE,UO,VO , DPOR ,P,P1,P2,P3,P4,PG1,PG2, 

5 MODE ,NOPNT , INPRO, KKONT ,MDI0G, KF,KFF, KPNT, SIG, ALPRM 
EOUIVAL ENCE (X,R) 

LF = 1 

IF (MT.NE.O) LF=2 
IF(K-3> 50,3,80 
50 IF(K-l) 2,1,2 

c K=1 PRINT INITIAL VALUES, HEAOINGS 

1 IF <MT) 13 , 10,11 

10 PRINT 200 
PRINT 206 

12 PRINT 201, XI, 0XI,0XM, XC,EMIN,EMAX 
DO 6 1=1, NT 

PRINT 208, I, Y ( I) ,1,7(1) 

6 CONTINUE 
PRINT 209 
LO=il+NT 
GO TO 60 

11 PRINT 202 
PRINT 206 
PRINT 206 
GO TO 12 

C K= 2 PRINT INTEGRATION STEP 

2 IF (LC + NT-58) 20,20,70 

70 GO TO (71,72) ,LF 

71 PRINT 200 
PRINT 209 
LC = 6 

GO TO 20 

72 PRINT 202 
LC = 6 

PRINT 209 

20 PRINT 203, X, OX ,P , OPOR , KKONT 
00 4 1=1, NT 

PRINT 205, 1, Y (I) , I , Z < I) , I,E< I ) 

4 CONTINUE 
PRINT 206 
LC-LC+2+NT 
GO TO 60 

C K= 3 PPINT FINAL CUTOFF INTEGRATION STEP 

3 IF (LC + NT-5fl)30 ,30 ,73 

73 GO TO ( 74,75) , LF 

74 PRINT 200 
GO TO 30 



75 °RINT 202 
*0 PRINT 211 

IF (MDIOG.NE.O) GO TO 59 
PRINT 204, XC 
00 5 1=1, NT 

PRINT 210, I, Y(I), I,Z(I) 

5 CONTINUE 
PRINT 207 
5 q CONTINUE 
60 RETURN 

K= 4 INOEX PAGE, PRINT HE AOINGS 

80 GO TO (81,«2) , LF 

81 PRINT 200 
PRINT 209 
LC = 6 

GO TO 60 

82 PRINT 202 
PRINT 209 
LC = 6 

GO TO 60 

200 FORMAT (56H1 RUNGE-KUTTA-GILL 7 AD A MS-MOULTON INTEGRATION ROUTIN 
IE) 

201 FORMAT (20H INITIAL C0NDITI0NS/16X1HX,17X5H0EL X,12X9HMAX DEL X, 

1 1 3 X , 5HCUT X,12X,9HMIN ERROR, 1 0 X , 9HMA X ERROR 75 X ,6E19. 7 ) 

202 FORMAT (40H1 RUNGE-KUTTA-GILL INTEGRATION ROUTINE) 

203 FORMAT ( ?E 1 8. 1 0 , 1 0 X , E 15 . 7 , 5X , E15 . 7 , 1 0 X ,6HKK0NT = I5 ) 

204 FORMAT (734H FINAL RESULT AT A CUT-OFF X OF,E15.7) 

205 FORMAT ( 1 Q X , 2 H Y ( , I 2 , 4H ) = , El 5 . 7 , 4X5H0ERV ( , 1 2 , 4H) = ,E15.7,4X, 

16HERR0R ( , 12 ,4H) = ,E15.7) 

206 FORMAT (1H ) 

207 FORMAT (1H171H1) 

208 FORMAT ( 9 X , 3H Y I ( , 1 2, 4H) = ,E15 . 7 , 3 X , 6H0ERV I ( , I 2 , 4H ) = ,E15.7> 

209 FORMAT ( 71 2X1 9H IN TEGR ATI ON RESUL TS/11 X , 1H X , 1 3X , 5HDEL X,23X,1HP, 

1 19X,4HD°DP/2 0X ,4HY Cl), 

22 3 X , 7H0ERV ( I) , 23X , 8HERROR f I ) / > 

210 FORMAT ( 10X,2HY(,I2,4H) = , E 1 5 . 7 , 4X5H0ERV ( , 12 , 4H) = ,El5.7) 

211 FORMAT ( 1H1) 

END 
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SUBROUTINE RKGI 


THIS ROUTINE USES THE RUNGE-KUTTA -GILL (R-K-G) METHOD OF 
INTEGRATION TO DETERMINE THE STARTING POINTS REQUIRED FOR 
USE WITH THE A-M METHOD 

COMMON QS<13,13) ,F (13,13) ,G< 13,13) ,H( 13,13) , AC (‘ 3) ,8C( 13) , 
Y(12),Y1(12),Y2(12I ,Y3(12) ,Y4(12) ,Z(12> ,Z1(12) , 
Z2(12) ,Z3(12> ,Z4<12) , YI(12) ,E(12) , X I , XC , X , XI , X 2 , 
X3,X4,0XI,DX,0XM,EMIN,EMAX,K0NT,LC,MP,MC,Mi,M2, 
M3,M4,MT,NT ,RE , UO , VO, DPDR , P , P 1 , P2 ,P3 , P4 , PG 1 ,PG 2 , 
MODE ,NOPNT , INPRO, KKONT , MOIOG, KF,KFF, KPNT ,SIG , ALPRM 
EQUIVALENCE <X,R) 

DIMENSION AA( 24) 

D1=SQRT (3.5000) 

02= SQRT (2.00 00) 

03=3. 0*01 
C 1=1 .0-01 
02 = 2 . 0-02 
C3=2. 0-03 
C4=1.0+Ol 
C5=2.0*02 
06=2. 0^03 
00 10 1=1, NT 
A A ( I ) =Z ( I ) 

CONTINUE 

DO 1 1=1, NT 

Y(I)=Y(I) +OX*Z(I) /2.0 

CONTINUE 

X= X1 + 0X/2. 0 

CALL OERV 

DO 2 1=1, NT 

Y ( I) =Y( I) +C1*DX*(Z(I) -AA (I) ) 

A A (I) =C2*Z (I) -C3* A A ( I ) 

CONTINUE 
CALL OERV 
DO 3 1=1, NT 

Y(I)=Y(I) +C4*DX*(Z(I) -AA(I)) 

AACI)=C5*7 (I) -C6* A A ( I ) 

CONTINUE 
X= X14-OX 
CALL OERV 
00 4 1=1, NT 

Y ( I ) = Y ( I ) +0X*Z(I) /6 • 0 -0X*A A ( I ) /3.0 

CONTINUE 

CALL OERV 

P =P1+SI G *OPDR*OX 

RETURN 

END 




SUBROUTINE ADMI 

r this ROUTINE UTILIZES THE ADAMS-MOULTON (A-M) preoictor- 
c CORRECTOR METHOD OF INTEGRATION OVER A FIXEO INCREMENT OF THE 
r INDEPENDENT VARIABLE IN ORDER TO EVALUATE THE DEPENDENT 

c VARIABLES Y(I) AT THE VALUE OF THE INDEPENDENT VARIABLE X*OX 

C 

COMMON QS (13, 13) ,F (13,13) ,G(13,13) ,H ( 13,13) , AC (13) ,BC(13> , 

1 Y(12) , Yl(12) ,Y2(12> ,Y3(12) , Y4 ( 12) , Z ( 12 ) ,Z1(12) , 

2 Z2(12) ,Z3(12) ,Z4<12) ,YI(12) ,E(12> ,XI,XC,X,X1,X2, 

3 X3,X4,0XI,DX,0XM,EMIN,EMAX,K0NT,LC,MP,MC,M1,M2, 

4 M3,M4,MT,NT,RE,U0,V0,DP0R,P,Pl ,P2,P3,P4 ,PG1,PG2 , 

5 MODEtNOPNT , INPRO, KKONT , MDI OG , KF , KFF , KPNT , S IG , A LPRM 
EQUIVALENCE (X,R> 

DIMENSION YY ( 24) 

IF(KF) 50,51,50 

50 IF (KFF) 52, 53, 52 

52 <F=0 
KFF=0 

GO TO 51 

53 KFF=1 

51 CONTINUE 

C OBTAIN TENATIVE VALUE OF NEXT POINT VIA AOAMS-MOUL TON INTEGRATION 

C CHECK ERROR FOR POSSIBLE CHANGE IN STEP-SIZE 

X=Xl+OX 
NM= 0 

Dl=0X/24. 0 
02=19.0X270.0 
00 5 1=1, NT 

Y(I)--Y1(I) +01*(55.0*Z1(I)-59.0*Z2(I)*37.Q*Z3(I)-9.0*Z4(I) ) 

YY ( I) = Y ( I ) 

5 CONTINUE 
CALL DERV 
00 6 1=1, NT 

Y(I) = Y1(I) +D1* <9. 0*Z(IM-19.0*Z1(I)-5.0*Z2(I>*Z3(I)) 

6 CONTINUE 

Call oerv 

C PRESSURE CALC FOR EACH REGULAR A-M STEP 

QQ = X2M X1*X1-X*X) +X1*(X*X-X2*X2) +X*(X2*X2-X1*X1) 

IF (QQ.EQ. 0.0) QQ=1.0E-16 

S1=(PG2*(X-X1) + PG1* ( X 2-X ) +DPDR* ( X 1-X2 ) )/QQ 
S2- (OPDR-PGl) / (X-Xl) -S1MX + X1) 

S3=PG2-S1*X2*X2-S2*X2 

p=P1 + SIG * (S3* (X- XI) +S2/2.*(X*X-X1*X1) +S1/3.* ( X*X*X-X 1 * X 1*X1) ) 

00 7 1=1, NT 

E ( I ) = ABS (02* (Y(I)-YY (I) )) 

IF (EMAX-E ( I ) ) 45,1,1 
1 IF (E(I) -EMIN) 8,7,7 

45 MM=MM+1 
GO TO 7 

8 M3=M3+1 

7 CONTINUE 

IF (MM) 46,46,3 

46 IF (M3-NT ) 10,2,2 

c oduSTlE Integration step 



2 IF <MC) 10,21,10 

21 IF CSIG *<X+4.*DX-XC) > 22,22,10 

22 IF (SIG * ( DX-* DX-DXN) ) 23,23,10 

23 CONTINUE 

IF (LC+NT-5R) 230,230,231 

231 CALL PRNT (4) 

230 IF (MP.EQ.O) GO TO 233 

IF (MOIOG.NE.O) GO TO 232 
233 PRINT 201 

CALL PRNT ( 2) 

PRINT 204 
LC=LC+2 

232 CONTINUE 
OX=DX+DX 
X1 = X 

P1 = P 

PG1-DPOR 
00 25 1*1, NT 
Z3 < I ) =Z4 < I ) 

Z1(I)=Z(I) 

Y1 (I) =Y (I) 

25 CONTINUE 
KK0NT=KK0NT+1 
CALL PKGI 

DO 26 1*1, NT 
E < I> =0 . 0 

26 CONTINUE 

IF (MOIOG.EQ.O) CALL PRNT(2) 

M 1 = 1 

GO TO 13 

C HALVE INTEGRATION STEP 

3 IF (M2) 30,40,30 

30 CONTINUE 
M2=l 

31 IF (LC + NT-5 8) 350,350,351 

351 CALL PRNT (4) 

350 IF (MP.EQ.O) GO TO 353 

IF (MOIOG.NE.O) GO TO 352 
353 KFFP1=KFF+ 1 

GO TO (451,452), KFFP1 

451 PRINT 202 
GO TO 453 

452 PRINT 206 

453 CONTINUE 
CALL PRNT ( 2 ) 

PRINT 204 
LC*LC+2 

352 CONTINUE 

IF (Ml) 34,32,34 

32 IF (SIG MX+OX+DX-XC) > 33,33,10 
33 IF (KFF ) 36,35,36 

35 OX*OX/2.0 

00 42 1*1, NT 
Y ( I) *Y1 (I) 

z<T)=zi(ir 



42 CONTINUE 

KKONT=KKONT+1 
CALL RKGI 
00 38 1=1, NT 
E ( I) = 0 « 0 
CONTINUE 

IF (MOIOG.EQ.O) CALL PRNT(2> 

K0NT=2 
KF = 1 
KFF = 0 
GO TO 13 
34 X1=X2 

P1 = P2 
M 1 = 0 

00 39 1=1, NT 

yi(i>=rzm 

Zl (I)=Z2 (I) 

39 CONTINUE 
GO TO 35 

10 H 2= 1 

IF (Ml) 11,13,11 

11 *1 = 0 
13 M3 = 0 

RETURN 
36 KFF=Q 
X 1 = X4 
Pi=P4 

DO 41 1=1, NT 
Y < I ) = Y4 ( I ) 

Y 1 <I>=Y4(T> 

Z ( I ) = Z4 < I ) 

Zl CI> =Z4 (I) 

41 CONTINUE 
GO TO 35 
40 Z2 ( 1) =1 • E12 
GO TO 13 

201 FORMAf ( 51H DOUBLING THE INTERVAL, ERROR ARRAY AS FOLLOWS) 

202 FORMAT (50H HALVING THE INTERVAL, ERROR ARRAY AS FOLLOWS) 

203 FORMAT < 8X ,6H ERROR < * I 2»5H) = ,E15.9) 

204 FORMAT (1H ) 

205 FORMAT ( 1 HI ) 

206 FORMAT (53H RE-HALVING THE INTERVAL, ERROR ARRAY AS FOLLOWS) 
END 
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SUBROUTINE DERV 

5 THIS ROUTINE EVALUATES THE DERIVATIVE VALUES AT EACH INCREMENT 
c OF THE INDEPENDENT VARIABLE 

C 

COMMON QS(13,13),F(13,13>,G(13,13> ,H(13,13) ,AC(13) ,BC(13) , 

1 Y(12), Yl(12> ,Y2(12) ,Y3(12) , Y4 < 12) , Z (12 ) ,Z1(12) , 

2 Z2(12) ,Z3(12> , Z4(12),YI(12> ,E(12I ,XI,XG,X, XI, X2, 

3 X3,X4,DXI,0X,DXM,EMIN,EMAX,K0NT,LC,MP,MC,M1,M2, 

4 H3,M4,MT,NT,RE,U0, V0,DPDR,P,P1,P2,P3,P4,PG1,PG2, 

5 MODE tNOPNT , INPRO, KKONT , MOI OG , KF ,KFF , KPNT , SIG 
EQUIVALENCE (X,R> 

T1 1-1. 0/2.0 
T1 2=1. 0/(3. 0*8.0) 

T13=l. 0/(5. 0*32.0) 

T14=l. 0/(7. 0*128.0) 

T15=1. 0/(9. 0*512. 0) 

BET Al = -UO/ (2. 0*R*R) 

MODD-MODE/ 2-1 
GO TO (100 ,200,300) , MOOD 
C 

c THIS SECTION IS FOR FORTH OROER POLYNOMIALS FOR U ANO V 
C 

100 CONTINUE 

IF (KONT.EQ.l) GO TO 120 
AC (1) =Y( 1) 

AC ( 3) =Y ( 2) 

AC(5)=-16.0*Y(1)-4.0*Y(2) 

BC ( 1) = Y ( 3 ) 

BC (3)=Y ( 4) 

8C(5)=16.0*(R-Y(3) )-4.0*Y(4) 

120 CONTINUE 
BET A3=0 . 0 
BET A4 = 0 • 0 
DO 140 1 = 1, 5, 2 
DO 140 J=1 ,5, 2 

8ETA3=BETA3+H(I,J) *BC(I) *BC(J) -G ( I , J) * AC < I ) * AC ( J) 

BETA4 = BETA44-F (I,J)*AC (I) *BC( J) 

140 CONTINUE 

BET A2 = - Y ( 1) *Y(3)/R + 2. 0*Y(4)/RE 

BETA3=BETA3/R+6.0*UO/ (R*RE) Ml. 0*Y(2) fAC(5) /2.0) /RE 
8ETA4 = -2»Q/R*BETA4M1,0*Y(4) ♦BC ( 51 /2 . 0 ) /RE 
T 24= Y ( 1 ) 

T31=Y(l)*(G(l,l)+H(i,l))^Y(2)*(G(l,3>+H(3,l>)+AC<5)*(G(l,5)+H(5,l) 

1) 

T32 = Y(1) MG(3»1)+H(1,3) ) +Y<2) *(G(3,3) ♦HO, 3) > ♦ AC (5) * ( G (3 , 5 ) +H (5 , 3 ) 
1) 

T 33=Y ( 1 ) *(G(5,1)+H(1,5) )+Y(2) *(G(5,3) ♦HO, 5) ) ♦ AC ( 5) * ( G ( 5 * 5) +H (5 , 5 ) 
1) 

T41=F (1, 1) *Y( 3>+F (1,3) *Y(4)*F(i,5) *BC(5) 

T42=F(3,1)*Y(3)+F ( 3 , 3 > * Y ( 4 ) +F < 3 , 5 ) *8C (5) 

T43=F(5,1) *Y(3)+F (5,3) *Y ( 4 ) +F (5 , 5 ) *BC (5) 

T44=F (1, 1) *Y( 1 ) + F ( 3,1) *Y(2)+F (5,1) *AC (5) 
T45=F(l,3)*Y(l)+F(3,3)*Y(2)+F(5,3)*AC(5) 

T 46=F (1,5) *Y(1)+F (3,5) *Y(2)+F (5,5) *AC(5) 



7(1) = (BETA 1*(T32-4.0*T33)-BETA3MT12-4.0*T13) ) /( (Til- 16. 0*T13) 
1MT32-4. 0*T33)-(T12-4.0*T13)* (T31-16. 0*T33) ) 

Z ( 2) = (BETA 1-(T 11-16. 0*T13) *Z ( 1) > / (T12-4. Q*T13> 

Z (3)=8ETA2/T24 

Z(4)=((BETA4-16.O*T46)-(T41-16.0*T43) *Z ( 1) - ( T42-4 . 0*T 43) *Z(2) 
1-(T44-16.0*T46)*ZC3) ) / (T45-4 . 0*T46 ) 

DPOR=Y( 3) **2/R+2. 0*Y(2) /RE-Y ( 1) *Z ( 1) 

RETURN 

THIS SECTION IS FOR SIXTH ORDER POLYNOMIALS FOR U AND V 
200 CONTINUE 

IF (KONT.EQ.l) GO TO 220 
AC (1)=Y (1) 

AC(3)=Y (2) 

AC (5) = Y (3) 

AC (7)=-64. 0*Y(1)-16.0*Y(2) -4. 0*Y(3) 

8C ( 1) =Y ( 4) 

8C (3) =Y ( 5) 

BC (5) = Y (6) 

BC(7)=64.0*(R-Y(4) )-16.0*Y(5) -4.0*Y(6) 

220 CONTINUE 
BET A5=0 • 0 
BETA6=0.0 
DO 240 1=1,7, 2 
DO 240 J=1 ,7, 2 

8ETA5=8ETA5>H ( I , J ) *BC (I ) *BC ( J ) -G ( I , J) *AC ( I ) * AC ( J) 

BET A6=BETA6+F (I , J) *AC (I)*BC(J) 

24C CONTINUE 

BETA5 = 8ETA5/R*-6.0*UO/ (R*RE ) ♦ ( AC ( 3 ) ♦AC (5) /2. 0^3 . 0*AC (7) /16.0) /RE 
BETA6 = -2.0*8ETA6/RMBC(3)^BC(5)/2. 0 + 3.0*BC(7)/16.0)/RE 
8ETA2=(-Y(1)*Y(4)/RV2.0*Y(5)/RE> /Y(l) 

3ETA3=12.0»Y(3) /RE + 2 . 0/R* ( Y ( 1 ) * Y ( 2 ) ♦ Y <4)*Y(5) ) 

BETA4=12. 0*Y(6) /RE+(Y (1) *Y (5 ) -Y (2 > *Y (4) ) /R 
T25=l • 0 
T3l=-Y (2) 

T32 = Y <1_) 

T41=-2. Q*Y (5) 

T45=Y ( 2 ) 

T46=Y(1) 

T51=Y(1) MG(1, 1)>H(1, 1) ) +YC2) ♦ ( G ( 1 ,3 ) +H <3 , 1 ) ) ♦Y (3) MG (1,5)>H(5,1) ) 
1*AC(7)MG(1,7I*H(7,1) ) 

T52=Y (1) MG(3 ,1)+H(1,3) ) + Y(2) * ( G ( 3 , 3 ) +H ( 3 , 3 ) ) ♦ Y ( 3) M G ( 3 , 5) +H ( 5 , 3) ) 

1 ♦AC ( 7 ) * ( G ( 3 , 7 ) +H ( 7 , 3 ) ) 

T53 = Y(i) MG(5,1)^H(1,5))M(2) MG(5,3)^H(3,5))M(3)MG(5,5)+H(5,5)) 
1^AC(7)MG(5,7)+H(7,5) ) 

T54=Y (1) MG(7,1)+H(1,7) )+Y(2> MG(7,3)+H(3,7))+Y(3)MG(7,5)+H(5,7>) 
1+ AC ( 7 ) * ( G ( 7 ,7 ) +H ( 7,7) ) 

T61=F(1,1) *Y(4)+F(1,3)*Y(5)+F(1,5)*Y(6)^F(1,7) *BC(7) 

T62 = F (3,1) *Y( 4 ) ♦F (3,3) *Y (5 ) +F (3 » 5 > * Y ( 6 1 ♦F (3,7) *BC (7) 

T63=F (5, 1) *Y(4) +F ( 5 , 3 ) *Y (5 ) +F ( 5 , 5 ) *Y ( 6) ♦F ( 5 , 7 ) *8C(7) 

T 64=F (7,1) *Y(4)+F (7,3)’Y(5)+F(7,5) *Y (6) +F ( 7 , 7 ) *8C ( 7) 

T65=F(1,1) *Y( 1)^F (3,1) *Y( 2) +F (5,1) »Y( 3) +F (7,1) *AC (7) 

T66=F(1, 7) »Y(1)+F (3,3)*Y(2)+F (5,3) ♦Y(3)+F (7,3) *AC (7) 

T67=F (1,5) *Y( 1) ♦F ( 3, 5 ) *Y ( 2) +F ( 5 , 5 ) *Y(3)^F(7,5) *AC(7) 

T68=F(1,7) *Y(1)+F(3,7)*Y(2)+F(5,7)»Y(3)+F(7,7) ♦AC(7) 



GAM1=T11-64.Q*T14 

GAM2=T12-16.0*T14 

GAM3=T13-4.0*Ti4 

ALPHl=T51-64. 0*T54 

ALPH2=T 52-16. 0*T54 

ALPH3=T53-4.0*T54 

FT1=0ETA1*ALPH3-8ETA5*GAM3 

FT2=G ANZ* ALPH 3-ALPH2*GAM3 

FT3=GAM1*ALPH3-ALPH1*GAM3 

Z(1>=<FT1*T32-FT2*8ETA3> /CFT3*T32-FT2*T31) 

Z(2> = <BETA3-T31*Z (1> ) /T32 
Z(3)=<8ETA1-GAM1*Z(1) -GAM2*Z<2> ) /GAM3 
Z(4)=BETA2/T25 

Z(5) = (BETA4-T41*Z C1>-T45*Z(4) ) /T46 

Z (6> = < <BETA6-64.0*T68 ) - (T61 -64. 0*T64> *Z ( 1 ) - ( T62-16 . 0 *T64> *Z(2>- 

l(T63-4. 0»T64) • Z (3 ) - C T 65-64 . 0*T68 > * Z ( 4) -{ T 66-16 . 0 *T68 ) *Z<5> ) / 

2 (T 67-4, 0 *T68) 

0P0R=-Y<1> *Z< 1>+Y(4>**2/R+2.0*Y(2)/RE 
RETUPN 

C 

c THIS SECTION IS FOR EIGHTH ORDER POLYNOMIALS FOR U ANO V 

C 

300 CONTINUE 

IF (KONT .ECU1> GO TO 320 
AC C 1) -Y ( 1) 

AC (3) -Y ( 2) 

AC(5)=Y(3) 

AC(7)=Y<4) 

AC(9>=-256.*AC<1> -64. * AC ( 3) -1 6. * AC (5 ) -4. * AC (7 ) 

BC ( 1> =Y < 5) 

9C (3) =Y ( 6 ) 

BC (5>=Y (7) 

BC ( 7) = Y ( 8 ) 

BC (9) =256. 0*(R-Y( 5> ) -64 . 0*Y ( 6 > -16 . 0* Y ( 7) -4. Q*Y (8) 

320 CONTINUE 
BET A7= 0 . 0 
8ETA8=0. 0 
OO 340 1=1,9, 2 
00 340 J=1 ,9, 2 

BETA7 = 8ETA7fH (I,J) *8C (I) *8C ( J) -G < I , J > *AC ( I > * AC < J) 

BET A8 = 8ET A8+F ( I , J ) *AC ( I ) *BC ( J ) 

340 CONTINUE 

BETA2=-Y(1) *Y(5)/R+2.0*Y<6>/RE 

BETA3=12. 0*Y< 3) /RE+2. 0/R* ( Y (1 ) *Y ( 2)*Y (5) *Y (6) ) 

8ETA4=12. 0*Y(7)/REMY(1)*Y (6 ) -Y < 2 ) *Y < 5) ) /R 

BETA5=2, 0*Y(5) * Y< 7 ) /R + 2. /3 . * Y < 2 ) * *2/R + Y (6 ) **2/R+4 . 0* Y ( 1) *Y (3) /R 
1+30 .0*Y(4) /RE 

BETA6 = -Y (3> *Y (5)/R-Y(2)*Y(6) /(3. 0*R) +3 . 0* Y ( 1 ) * Y ( 7 ) /R* 3 0 . 0*Y ( 8 ) /RE 
BETA7 = BETA7/R + 6.0*UO/ <R*RE > ♦ < AC ( 3 ) +AC <5) /2 . 0 + 3 . 0 * AC < 7 > /16 . 0 + 

1 AC < 9) /16.0)/RE 

BETA8 = -2 .0*8ETA8/R+(8C<3) + BC C5) / 2 . 0*3 . 0*BC < 7 ) / 16 . 0+8C ( 9) /16.) /RE 
T26=Y ( 1) 

T 3 1 = - Y < 2 ) 

T 32= Y ( 1 ) 

T41=-2. 0*Y (6) 



T46=Y (2) 

T47=r (1) 

T51=-3. 0*Y(3) 

T52=Y<2> 73.0 
T53=Y ( 1 ) 

T61=-4. 0*Y (7) 

T62=-2./3.*Y(6> 

T66=Y (3) 

T67=Y(2) 

T68=Y ( 1) 

T71=Y(1) *(G(1,1H-H (1,1) )+Y(2) *(G(i,3> +H(3,1> )+Y(3)* (G (1,5)+H(5,1) ) 
l + flC(7) MG (1,7) +H(7,1) ) + AC(9> * (G C 1 , 9) +H (9, 1) ) 

T72=Y(1) * (G(3,1)*H(1,3) > + Y(2> MG (3,3) +H(3,3) ) +Y (3 ) M G ( 3 , 5 > +H ( 5, 3 > ) 
1*AC(7)* (G(3,7> ♦H(7,3) )+AC<9)* (G (3 , 9) +H <9, 3 ) ) 
T73=:Y(1)MG(5,1)*H(1,5))+Y(2)MG(5»3)+H(3, 5) ) ♦Y ( 3) * (G ( 5 , 5 ) +H C 5,5) ) 
lfAC(7)MG(5,7)+H(7,5> > + AC(9) * (G ( 5 , 9) + H (9, 5 ) ) 
T74=Y(1)MG(7,1>*H(1,7)>*Y(2>MG(7 # 3>+H<3#7))+Y(3>MG<7,5> + HI5,7> ) 
1*AC(7)MG{7,7>+H(7,7) ) +AC(9) * CG (7 ,9) 4-H (9, 7) ) 

T75=Y(1> ♦ (G(9,l)+H(l, 9) >*Y(2) MG(9,3>+H(3,9))M(3>MG(9,5>+H(5,9>> 
i*AC(7) * (G(9,7) +H(7,9> ) + AC(9)* (G ( 9 , 9) +H (9, 9) ) 
T81=Y(5>*F(i,l)+Y(6)*F(l,3)+Y(7)*F(l,5)+Y(8)*F(l,7)+BC<9>*F(l,9) 

T82=Y(5) *F (3, U+Y (6>*F (3,3>*Y (7>*F (3,5>Vf (8) *F (3 ,7)+BC (9>*F (3 ,9) 
T83=Y(5> *F (5,1) *Y (6) *F(5,3)+Y<7> *F(5,5)+Y(8> *F(5, 7)+BC (9) *F(5,9) 

T8 4=Y(5)*F(7, 1 > +Y ( 6) *F ( 7 , 3) +Y (7 ) *F (7 , 5 M*Y ( 8 ) *F ( 7, 7) +BC (9) *F ( 7 ,9 > 
T85=Y(5 )*F(9, 1) +Y ( 6) *F (9 , 3) +Y (7) *F (9> 5>+Y < 8> *F (9 , 71+BC ( 9) *F (9 , 9) 
T86=Y(1)*F (1,1)^(2)*F(3,1)+Y(3) *F (5, 1) +Y (4) *F ( 7, 1) ♦ AC < 9) *F (9 , 1) 

T 87=Y ( 1 > *F ( 1, 3) ♦*? ( 2) *F < 3 , 3) +Y (3) *F (5 , 3) +Y ( 4) *F (7 , 3) ♦ AC (9) *F (9 , 3 ) 

T 8 8 = Y ( 1 ) *F (1,5) + v (2) *F(3,5)*Y<3) *F«5,5)+Y(4> *F (7,5> + AC (9) *F(9 ,5) 

T8 9=Y(1I *F (1, 7M-V ( 2> *F ( 3 , 7) ♦ Y ( 3) *F (5 , 7) +Y (4) *F (7,7) ♦AC (9)*F (9,7) 

T810=Y(1)*F(1,9)*Y(2) *F (3 , 9) ♦ Y ( 3) *F (5 , 9) + Y ( 4 ) *F ( 7 , 9) +AC (9) *F (9, 9) 

ALPHA 1=T11-256.0*T15 

ALPHA2=T12-64.0*T15 

ALPHA3=T13-16.0*T15 

ALPHA4=T14-4. 0 # T15 

GAMl=T71-256. 0*T75 

G A M2 = T 72-64 .0 *T75 

G AN3=T73-16 .0 *T75 

5AM4=T74-4.0*T75 

RH 01= ALPHA l*GAM4-ALt»HA4*GAMl 

RH02=ALPHA2»GAM4-ALPHA4*GAM2 

RH03= ALPHA 3*G A M4- ALPHA 4*GAM3 

EP l=RHOjL * T 5 3- R H03 * T5 1 

FP2=RH0 2*T53-P.H03*T52 

Z (1) = ( ( (8ETA1*GAH4-ALPHA4*8ETA7) *T53-RH03*BETA5) »T32- EP2*8ETA 3) / 
1(EP1*T32-EP2*T31) 

Z(2)=(8ETA3-T31*Z (1) ) /T32 

7(3) =(8ETA5-T51*Z(i)-T52*Z(2> ) 7T53 

Z (4) = (BETA 1-ALPHA 1*Z( 1) - ALPH A2*Z ( 2) - ALPHA3 * Z ( 3 ) ) /ALPHA 4 
Z (5)=8ETA2/T26 

Z(6)=(8ETA4-T41*Z (1)-T46*Z(5) ) /T47 

Z(7)=(9ETA6-T61*Z (1)-T62*Z(2) -T66*Z C 5 ) -T67* Z ( 6 ) ) /T68 
Z (8) = ( (BETA8-256. 0*T81O> -< T8 1 -256 . 0»T85 ) * Z ( 1 > - ( T8 2-64 . 0 *T85> * 1 (2) 

1- (T83-16«0*T85) * 1 ( 3) - (T84-4. 0*T85 ) *Z ( 4)- ( T86-256 . 0*T8 1 0 ) *Z ( 5) 

2- (T 87-64.0 *T810) * Z (6 ) - ( T88-16 . 0*T810 ) *Z ( 7 ) ) / (T89-4. 0 *T810) 

OPDR= Y ( 5 ) **2/R+2.0*Y(2)/RE-Y(l>*Z(l) 


RETURN 

END 


SUBROUTINE VECTOR ( VEL MAG , VEL ANG) 

COMMON QS(13, 13), F (13, 13) ,G(13,13> ,H(13,13> ,AC(13) ,BC (13) , 

1 Y(12),Y1(12> ,Y2(12J ,Y 3 ( 12) , Y4 ( 12) , Z ( 12) , Z 1 ( 12 > , 

2 Z2(12) ,Z3(12> , Z4(12),YI<12) ,E(12) , XI , XC , X , XI , X2 , 

3 X3, X4,0XI,DX,DXM,EMIN,EMAX ,KONT,LC,MP,MC,Mi,M2, 

4 M3 , M4, MT , N T , RE , UO , VO , DPDR , P , PI , P2 , P3 , P4 , PG 1 , PG 2 , 

5 MOOE,NOPNT,INPRO,KKONT,MDIOG,KF,KFF,KPNT,SIG,ALPRM, 

6 ATES<64,65> ,XTES(64) , 

7 RSV (31) ,ETASV(31) ,PTSV(31) ,PSV(31) ,TQUESV(31) , 

8 ORORSV (31) , AC0(13) ,ZINF(13) ,ZVUV(21,13) ,ZVW<21, 13) , 

9 I NUM (31) ,RINF(11> ,AUP(7,3) ,UH(21) ,VW(21) ,RI1, 

1 ZW(21) ,HWC21) ,Q(13,13> ,XPNT(31) ,RI2 , 1 S , LU , NTP , P I , NTH , 

2 VM AG (3 1 ) , VANG (31) 

EQUIVALENCE <X,R) 

DELZ=0.5/FLOAT (LU-1) 

LIM=(LU-1) /2 

UW ( 1) = AC ( 1 ) 

V W ( 1) =8C ( 1 ) 

WW(1)=0.0 

00 521 1-2, LU 

UW(I)=0. 0 

VW ( I > = 0 • 0 

WW ( I ) =0 • 0 

00 518 J-l ,NTP , 2 

UW(I)=UW(I)+AC(J) *ZVUV(I,J> 

VW(I)=VW(I)+BC(J) *ZVUV(I,J) 

FJ = J 

WW (I)=WW (I )-( AC(J) /R*ACO(J> > *ZVW(I ,J) /FJ 
518 CONTINUE 
521 CONTINUE 
UA VG= 0 .0 
VAVG=0.0 
DO 20 1=1, LIM 
JJ=2*I 
JJJ=JJ+1 

UAVG=UAVG+4.*UW ( JJ)*2.*UW( JJJ) 

VA VG=VAVG+4.*VW(JJ)+2.*VW( JJJ) 

20 CONTINUE 

UAVG= (UAVG>UW(1)-UW(LU) ) *OELZ/3. 

UAVG=UAVG/0.5 

VAVG= (VAVG+VW (l)-VW(LU) ) *OELZ/3. 

VAVG=VAVG/0.5_ 

VELMAG=SQRT(UAVG*UAVG+VAVG*VAVG) 

ARG=-VAVG/UAVG 

VELANG=ATAN(ARG)* 180. /3. 1415926 
RETURN 
E N 0 
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BOYACK OUTPUT DATA SHEETS 
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SUMMAPY 

INCOMPRFSSIPLF FL )H '>c T rF<-N ? 0 T A T T NG DISKS 

POYACK INTEGRAL M:T-»'V* v * A r C H ■” D TO POYH F-0 START-UP 

PTH 0 °OER SOLUTION tJMTF ORH INLET ORQFILFS 


un - 

- . * n 2 0 U 

' ,n = 1.13 3 0 

RF = 

1 . 0 0 0 0 


p 

ETA 

n t 

P 

9P9P 

T 

1.-30 

0 . C C 0 ] 

. c '• *117 

0.000000 

1 .086893 

0.000000 

. Q 5 0 

• 26 7 0 

• 6 c 8 1 6 4 

. 062194 

1.213264 

.022051 

.HOC 

. 3 7 76 

. ^ l r«05 

.121975 

1.177997 

.033653 

. 85 r 

. 45«4 

. 7 7 “ 7 2 

.180002 

1.144453 

. 044638 

. «0u 

. 52 54 

. c <**76 7 

.236397 

1.113023 

.054987 

. 75 C 

.5812 

.-•71^5 

.291325 

1.084053 

. 064716 

. 70 G 

.5252 

. j 4 i • 7 4 R 

. 344879 

1.058166 

.077818 

. 650 

. f 502 

, p o 7 n 9 9 

.397228 

1 . 0 36C 09 

. C 8 2 2 9 1 

. 60 C 

.6*76 

1. 44729 

. 448459 

1.018645 

.090119 

.55 0 

. 7083 

1 . * R4R2 0 

. 49 Q 0 5 0 

1.007319 

.097333 

. 5U G 

. 7271 

t . 1 452R8 

.549427 

1.003891 

.103941 

. 45 C 

. 7724 

1.1 ° f; 5:ll 

.599641 

1.011174 

.109890 

• 4 C C 

. 73 65 

1 . •»4<=-5?i 

.650651 

1.0 33352 

.115223 

. S5i 

. 73 55 

1 . 2 9 R 9 - 3 4 

. 703364 

1 .077522 

.119936 

. 30 0 

. 72 R 7 

1 . tir 4 8,26 

. 758956 

1.154629 

.124014 

. 25 0 

. 71 75 

1 . 415547 

.819677 

1 . 2 8 7 7 1 0 

.127465 

. 200 

. 69«0 

1 «4®5 47 4 

. 889604 

1.526034 

.130298 

. 15 ; 

.6717 

1 . 5 7 1720 

.975850 

1.996221 

.132489 

.100 

.6206 

1 * 6 °6 62 3 

1.100750 

3.209949 

.134066 

. 05 0 

. 5478 

1 . or, 7729 

1 . 767858 

9.723533 

.135015 


NO U 

prof I L 1 T NFL r - O r r pm 



TIMF 

IN EXECUTION = SECONDS 

COC 6400 

12 720/71 


INTEGRATION S TF pc i,oy 
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SUMMARY 

INCOMPRESSIBLE FLOW BETWEEN ROTATING DISKS 

ROYACK INTEGRAL METHOD MATCHED TO BOYD F-D START-UP 

8 TH ORDER SOLUTION UNIFORM INLET PROFILES 


uo = 

-. C 200 

VO = 1.1010 

RE - 

2.0000 


0 

ETA 

PT 

P 

QPDR 

T 

1 . GOO 

C . 0 0 0 0 

.595113 

0.000000 

1.086893 

0.000000 

.050 

. 2619 

.652860 

. 056927 

1.097380 

.021453 

. 900 

• 

CD 

. 706651 

. 110718 

1.055673 

.033056 

. 85 0 

.4627 

.758383 

. 162450 

1.014972 

.044031 

.800 

. 5363 

.808189 

.212257 

.975396 

.054396 

. 750 

.5969 

. 855973 

. 260040 

.937312 

.064116 

. 700 

.6469 

.901969 

. 306037 

.900916 

.073224 

.650 

.6880 

.946125 

.350192 

.866745 

.081690 

.60 0 

.7216 

.9886 39 

. 3927Q6 

.835306 

.089530 

. 550 

. 7487 

1.029732 

.433799 

.807341 

.096754 

.50 0 

.7700 

1.06 947 2 

.473539 

.784068 

.103339 

.450 

.7859 

1.108198 

. 512265 

.767105 

.109297 

.400 

. 7968 

1.146161 

.550229 

. 759050 

.114612 

. 35 0 

. 80 30 

1.184276 

.588343 

.764147 

.119333 

. 30 C 

. 8041 

1.222868 

. 626936 

.789861 

.123400 

. 25 P 

. 79R9 

1.263813 ~ 

.667880 

. 851448 

.126864 

.20 0 

.7893 

1. 309389 

.713456 

.983178 

.129699 

. 150 

. 7 7 0 1 

1.364685 

. 768753 

1.280558 

.131891 

. 100 

. 7357 

1.445704 

. 849771 

2.152126 

.133465 

. 05 U 

.6515 

1 . 6 44169 

1 .048237 

7.871314 

.134433 


NO U PROFILE INFLECTION 

TIME IN EXECUTION = 73.66 SFCONDS COC 64QG 12/17/71 


INTEGRATION STEPS 


5465 



55 

SUMMARY 

INCOMPRESSIBLE FLOW BETWEEN ROTATING OISKS 

BOYACK INTEGRAL METHOD MATCHED TO BOYD F-D START-UP 

8TH ORDER SOLUTION UNIFORM INLET PROFILES 


UO = -.0200 


VO = 1.1050 


RE = 3.0000 


R 

ETA 

PT 

P 

DPOR 

T 

1.000 

0.0000 

.595113 

0 .000000 

1,086893 

O.OOOuuU 

.950 

.2547 

.651597 

.055572 

1.065974 

.020826 

. 90 C 

.3678 

.703931 

.107906 

1,022074 

• G32482 

.850 

.4590 

.753695 

.157670 

.979636 

.043403 

.800 

.5349 

.801909 

. 205884 

.937674 

.053818 

.750 

.5973 

.847746 

• 251 72(F 

.897115 

.063537 

. 700 

.6492 

. 891589 

.295564 

.857920 

.072629 

.650 

.6924 

.933517 

. 337492 

.820502 

.081095 

.601) 

.7230 

.973627 

.377602 

.785277 

.088934 

.550 

.7571 

1.011923 

. 415898 

.753001 

.096124 

.500 

.7806 

1.049056 

.453031 

.724253 

.102753 

• 45C 

.7988 

1.084524 

.488499 

.700955 

. 1 0 8692 

. 400 

.8120 

1.119115 

.523090 

.684968 

.114025 

. 350 

.8206 

1.153276 

.557250 

.680074 

.118746 

. 30 C 

. 8243 

1.187375 

. 591350 

.692928 

.122812 

. 250 

.8229 

1.222810 

.626785 

.736824 

. 126257 

.20 0 

.8153 

1.262062 

.666037 

.844165 

.129097 

.151 

.7990 

1.309625 

.713600 

1.107415 

.131297 

. 100 

.7662 

1.382341 

.786316 

1.961255 

.132898 

.050 

.6780 

1.573637 

. 977612 

8.036757 

.133869 


NO U PROF 

ILE INFLECTION 



TIMF TN 

EXECUTION = 58.52 SECONDS 

CDC 640 U 

12/21/71 


INTEGRATION STE°S 


3822 



31JMMAP Y 

INCOMPPESSIDLS FLOW 9ETWEEN ROTATING OTSKS 

99YACK TNTFGPAL METHOD MATCHED TO ROYD E-U START-UP 

8TH ORDER SOLUTION UNIFORM INLET PROFILES 


1)0 = 

-.0200 

v d = l.iono 

RE = 

4.0000 


p 

FTA 

PT 

o 

DPDP 

T 

1.000 

0 . u 0 0 0 

.595113 

U . OOU OOU 

1.086893 

U. 000000 

. 05 0 

. 2475 

.661553 

.055417 

1.056012 

.020232 

. OOu 

.7616 

. 7 03170 

.107034 

1.01072? 

.031901 

. 850 

.4 541 

.757391 

. 156455 

.967410 

.042873 

.800 

. 5304 

. 7 99869 

. 203733 

.924982 

.053219 

. 750 

. 59^9 

. 845198 

.249062 

.883472 

.062969 

. 700 

.6465 

.887200 

. 292064 

.843520 

.072031 

.65 0 

.6934 

. 929383 

. 333247 

.805057 

. 080498 

.600 

.7270 

.963819 

. 372683 

.768527 

.088362 

. 55b 

.7569 

1 . 0 06247 

. 4 1U11 t 

.734898 

.095551 

.500 

.7817 

1.042315 

. 446179 

.704625 

.102159 

.45 0 

.8004 

1.076982 

.480846 

.679309 

.108136 

.400 

.8145 

1.113330 

.514194 

.661227 

.113451 

. 350 

. 8240 

1.143737 

.547101 

.653763 

.118172 

. 301 

. 8786 

1. 175739 

.579602 

.663608 

.122215 

. 250 

. 8281 

1.210095 

.613959 

.705215 

.125706 

.230 

.821 3 

1.247461 

.651326 

.810793 

.128529 

. 151 

.80 57 

1.293495 

.697359 

1.081618 

.130734 

. 100 

. 7723 

1.365987 

. 769351 

2.002421 

.132345 

. U50 

.6761 

1.571918 

. 975782 

9.303497 

.133333 


U PROFILE INFLECTS AT P= .1)60 758 


TIME IN EXECUTION = ST. 20 SECONDS COD 6400 12/15/71 

TNfFOPATTOK OTEPC = 2035' 
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SUMMARY 

I NCOM P°E S S I P L E FLOW RETWEEN ROTATING DISKS 

flOYACK INTEGRAL METHOD MATCHED TO RDYD F-D START-UP 

*5 TH OROER SOLUTION UNIFORM INLET PROFILES 


uo = 

-.0200 

VO = 1.1010 

RE = 

6. OOuQ 


p 

ETA 

PT 

P 

OPOR 

T 

1. JOG 

l . 0 0 0 0 

.595113 

0.000000 

1.086893 

0.000000 

. 950 

. 2319 

.652010 

. 055612 

1.058647 

.018957 

. 900 

. 3494 

.703807 

.107409 

1.009887 

.030835 

. 85 G 

.4427 

.753180 

. 156782 

.966403 

.041813 

.800 

.5196 

.800439 

.204011 

.924000 

.052155 

.75 0 

.5833 

.845400 

. 249003 

.882812 

.061839 

. 700 

.6368 

.888648 

. 292251 

.842704 

.070959 

.650 

.6811 

.929662 

. 333264 

.804515 

.079395 

.601 

.7179 

.968824 

. 372426 

.768450 

.087210 

.55 it 

.7484 

1.006497 

.410099 

.734960 

.094448 

. 500 

. 7730 

1.042580 

. 446183 

.705269 

.101056 

.45 0 

.7923 

1. 0 77301 

.480904 

.680892 

.107035 

.40 0 

.8066 

1.111764 

.514367 

.664354 

.112352 

. 351 

.8162 

1.143675 

.547277 

.659619 

.117048 

. 30 P 

. 8209 

1.177236 

.580838 

.674864 

.121174 

. 251 

. 8201 

1.211453 

.615055 

.725386 

.124583 

.201 

. 8128 

1.250300 

.653902 

.853080 

.127427 

. 150 

. 7949 

1.301251 

.704854 

1.201184' 

.129710 

. 10 0 

. 7567 

1.383956 

.787559 

2.405336 

. 131325 

. 050 

.6377 

1.655404 

1.059007 

12.614015 

.132363 


U PROf-lLF INFLECTS AT R= .99292b 
U PROFILE UNINFLECTS AT P= .9*9690 
U PROFILE INFLECTS AT R= .131660 


TIME IN EXECUTION = 42.23 SECONDS " CDC 6400 12/17/71 


INTEGRATION STEP^ 


2127 



SUMMARY 

INCOMPRESSIBLE FLOW BETWEEN ROTATING DISKS 

boyack Integral method matched to boyo f-o start-up 
sth order solution UNIFORM INLET PROFILES 


uo = 

-.0199 

VO = 1.1000 

RE = 

8.0000 


R 

ETA 

PT 

P 

DPDR 

T 

1.000 

0.0000 

.595113 

0.000000 

1 .086893 

0. 000000 

.95 0 

.2182 

.653068 

. 056398 

1.070965 

.017863 

.900 

. 3373 

.705062 

. 108392 

1.019969 

.029007 

.850 

.4315 

.755101 

. 158431 

.976555 

.040838 

.801 

.50 79 

.802393 

. 205723 

.935259 

.051079 

.750 

. 5721 

.848259 

.251589 

.894676 

.060828 

.70 0 

.6256 

.892127 

.295457 

.855606 

.069947 

.651 

.6697 

.933813 

.337143 

.818668 

.078383 

.600 

.7067 

.974230 

. 377560 

.783702 

.086297 

.551 

.7365 

1.012223 

.415553 

.752586 

.093438 

.500 

.7609 

1.049408 

.452818 

.725250 

.100089 

. 45C 

.7797 

1.085062 

.488392 

.704383 

.106030 

.401 

.7934 

1.119806 

.523135 

.692564 

.111352 

.350 

.8023 

1.154730 

.558068 

.694943 

.116113 

. 30(1 

.0059 

1.189880 

.593210 

.721353 

.120192 

. 250 

.60 39 

1.227333 

.630663 

.792828 

.123656 

.201 

. 7945 

1.270462 

.673782 

.963252 

.126512 

.151 

.7730 

1.327623 

. 730953 

1.409520 

.128766 

. 101 

.7286 

1.428221 

.831551 

3.021832 

.130432 

. 050 

.5815 

1.805062 

1.208392 

18.210124 

.131557 


U PROFILE INFLECTS AT P= .992926 
U PROFILE UNINFLECTS AT R = .970992 

U PROFILE INFLECTS AT R= .275327 


time In execution = 38.&o~ seconds cd(T 6400 "" 12 / 21/71 


INTEGRATION STEPS 


1746 
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SUMMARY 

INCOMPRESSIBLE FLOW BETWEEN ROTATING OISKS 

BOYACK INTEGRAL METHOO MATCHED TO BOYD F-D START-UP 

«TH ORDER SOLUTION UNIFORM INLET PROFILES 


UO = 

-.0 199 

VO = 1.1000 

RE = 

10.0000 


p 

ET A 

° T_ 

P 

OPQR 

T 

1. uoo 

0 . 0U 00 

.595113 

0 .000000 

1.086893 

0.000000 

.950 

. 2054 

.653644 

. 056767 

1.086212 

.016819 

.900 

.2260 

.706353 

.109475 

1.033821 

.028848 

.850 

.42 0 v 

. 757090 

. 160213 

.990776 

.039902 

.800 

.4976 

.805419 

.208541 

.950123 

.050214 

. 750 

. 5620 

. 852357 

.255479 

.910458 

.060022 

. 701 

.6145 

.896477 

.299600 

.873274 

.069020 

.65 0 

.6589 

. 939633 

. 342756 

.837535 

.077566 

. 6 U * 

.694/ 

.980262 

. 383385 

.805265 

.085327 

.55 0 

.7244 

1.019899 

.423012 

.776273 

.092566 

.501) 

. 7482 

1.058430 

• 461553 

,752252 


.451) 

. 7662 

1.095451 

.498574 

.735764 

.105164 

.40". 

. 7790 

1.131.420 

.534543 

.729956 

.110425 

. 351 

.7869 

1. 168446 

. 571569 

.741110 

.115202 

. 301 

. 7894 

1.206201 

.609324 

.781260 

.119301 

.251 

.7859 

1 .247178 

• 650 3C 1 

.876942 

.122790 

. 201 

. 7742 

1.295559 

. 698662 

1.097116 

.125678 

. 15 C 

. 7490 

1.364265 

. 767387 

1.699254 

.128028 

. 10 0 

.6949 

1.490424 

. 893547 

3.923711 

.129738 

. 05 0 

. 5237 

1.995181 

1 .398303 

25.3 4 642 u 

. 130915 


U PROFILF INFLECTS AT R = .992926 


TIME I N EXECUTION = 36.42 SECONDS COC 6400 


INTEGRATION S TE D S = 15C9 


12/17/71 
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PERFORMANCE MAPS 
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1. The performance maps in this section are plotted for 
V Q = 1.0 with parabolic inlet profiles at the indicated location. 






103 





















08 


R. = 0.7 



0 4 8 


N 


E 






























115 







116 












119 



































127 










129 


2. The performance maps in this section are plotted for 
V Q = 1.1 with parabolic inlet profiles at the indicated R. location. 
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The performance maps in this section are plotted for 
th parabolic inlet profiles at the indicated R.. location. 
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4. The composite maps in this section are plotted for 
V Q = 1.1 with parabolic inlet profiles at the indicated location. 
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5. The graphs in this section are plotted for the indicated 
parameters with parabolic inlet profiles. 
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V = 1.3 
o 

U = - 0.04 
o 



R: 








U = - 0.02 
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U = - 0.02 
o 

N RE = 6,0 










OD vj OOiik. CJ KJ 
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6. The performance maps in this section are plotted for 

1/ = 1.1 with uniform inlet profiles at the indicated R. location, 

o r i 
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APPENDIX D 


DESIGN PROGRAM 



OC^OOClOOOfCl 


TURBINE DESIGN STUDY PROGRAM 

THIS PROGRAM USES boyack results as input and calculates turbine 

PERFORMANCE, 

from boyack results, uo.nre* efficiency* and Ri are entered, 
ALSO* the head (FEET OF WATER) IS ENTERED, 

THE PROGRAM YIELDS Q. OMEGA* POWER* NO_, OF DISKS.ETC, 


DIMENSION FLu I D < 5 ) 

1 READ 2 , FLU I D 
READ 3. FNU.PMH 
IF(RHO) 10 00 * 1000 , 5? 

5 READ 10 , UO , VO , RE * R I * ETA » PT • T 
READ 3.R0 
CON* 1.0 
D* , 015 
RMIN=1, 

RM A X * 24, 

RQBAR = RO 
HEAD=0 . 

15 HEAD * HEAD ♦ 20. 

I F ( HE AD-22 0 . ) 25.25,1 

25 PTBAR = HEAD / 2.31 

ALPHA s ATAn2(U0* VQ)*36Q./6,2831S5 

PRINT 30 

PRINT 40 

PRINT 60, FLUID 

PRINT 65,F.\U*RH0 

PRINT 70 , UO , yo * RE 

PRINT 75,RI,cTA,PT*T 

PRINT' 12* ALPHA 

print 11 * d 

P r R I NT 80, PTBAR* HEAD 
print 90 
PRINT 95 

IF (RO) 26.29,27 

26 CONTINUE 
DR*i. 

RQBAR = RMlN - DR 
120 ROBAR = ROBAR + DR 

IF ( ROBAR . GE . 2 , ) DRs2, 

IF (ROBAR - RMAX) 130, 13<J* 15 
?T CONTINUE „ „„ 

130 OMEG s PTBAR « 32.2 * 144, * 144, / (PT«RHO*ROBAR»RQBAR) 

OMFG = SORT (QMEG) 

OMp.GA = OMEG * 60, / (6.2831853Q72) 

H = RE * F*'U / OMEG 

H = SQRT(H) , . 

Q = 6,283lQ53o/2 * UO « RQBAR * ROBAR * OMEG * H * 60, / 1728, 
TBAR s T * R H 0 * ROBAR * ROBAR * R 08 AR * ROBAR * OMEG * OMEG • H 
TEAR a TBAP / (32 . 2»1728 , *144 , ) 

RIBAR * RI * ROBAR 



t_ T W 


CON* . 0 

2 A COM*CO»i. 

FL s ROBAR <* CON 
200 FN = FL / (0*H> 

TOR 3 FN * T0AR 
QQ s FN * Q 

POW = TOR * OmEG / 550 ; 

PRINT 110 # RQBaR » R I BAR » FL » H i FN # TOR , POW, QQ » OHEGA 
IF<C0N.EQ,2, } GO TO 28 
IF (RO) 28,26,29 
GO TO 24 
28 CONTINUE 
GO TO 120 
CONTINUE 
PRINT 41 
PRINT Rl 
PRINT 96 

P'I2s6, '2831*53072 


DPTsPTBAR/?0 • 

TOnTso 

PTOF a PTB AR+ ( Op T , ) 
36 PTOF s PTQF-t)PT 


K0NT»K0NT*1 
OHEG « PTQF 


32.2 * £44, * 144, / ( PT*RWO*RCB AR*RQ V 


OMF3 ■ SORT {qMpG) 
■Q = 6V2B31S53ff7g • 


* UO • * RO BAR * TTmEG « H * 3UY / VW& s 

TB A T s T # RHQ •» ROBAR # ROBAR * ROBAP * ROBAR * OMEG * OMEG • 


TSAR a TBAR / (32.2*1728**1*4.) 

TOR 3 FN * T 6 a R 

FOW r 3 TOR * OMEG / 550*. 


OMEGA s OMEG * 60, / (6,2831853072) 
EFF57ffAR*0'MEG*6 0./<Q*PTOr*144, ) 


QQ=FN«Q 

PRINT ill, QQ »OMEGA, Pnw.foR*EFF,pTOF 
IF( KONT.EQ.il) GO TO 15 
GO TO 36 
1C00 CALL EXIT 

2 FORMAT’ (5a5) 

3 F0RMAT(2F1?.G) 

10 format ( 7Fio,0) 

11 F0RMAT(5X,4N d 3 R8,4#2X,2hIN) 

12 FORMAT < 5X , 15M NOZZLE ANGLE *F6.2»5H DEG.) 

‘40'TORMAT(5X^35M MULTIPLE DISK TURBINE DETTGN STUDY",/,?*, 2* w BASTE O’ 'W 
1 8QYACK RESULTS,//) 

60 F0RMAT(5X,5A5) 

65 F0RMAT(5X,22 h K 1 \£ k a T I C VISCOSITY »F9 . 4 , IX , 8H 1 NSQ/SEC , / » 5X , 10H DEI* 
IS X TY aF9.4,6H LB/FTCU,/) 

70 F 0 RMAT( 5 X, 5 W u° =F6.2,3X,gH VO sF6.2»3X#5H RE sF6.2) 

75 F0RMAT(5X.5H pi =F7.3,2X.6H~ETA «Fff,3','2Tr»1WT *F10 ,6 * 2X» 3KT sFlO.fc 


30^FORfiAT(5X,BM pTB AR =F 8 ", 2 . IX , 8HLBF/ INSQ , 3X , 6HWE AD »F8 . 0 , IX , 2KFT • / 5 

90 FORMAT ( 7K , 2HRC » 8X, 2HR I ,'5X, lHL*6X, 1HH. 5X» 1HN, 5X» 3WT0R»6X» 2HHP, 3X, 6h 



lVnt FL#5X,3 h r FM) 

V5 rORMAT(«X» 4H( IM) > 3X,4HeI f 'l) *3X»4H( IN) ,3X', 4>T(TN) »7y,7ff(TT-l.t)»8X71C^ 
KFTCU/MIN) »/> 

110 F0RM*TC 5X » F5. 2 » 2 x , F5 , 2* 2X # F5 . 2 * 2X # T5 . 4 »1X#T5. ff # IX » V7 , <i , £X* F7 , 2 ,TX . 
lF7,2*2X»^6.0) 

41 F0R^AT£//>5X»21H OFF DESIGN OPERATION »/» 5X.17HWJ TW HEAD VAffYTNGT,/) 
91 F0RHAT(5X,6HV0L F L • 5X . 3HRRM i 6X , 2HHP » 6X . 3HT0R » 5X . 3HJET kj 3X»5H 

1PTBAR )' ' ~ 

96 FORMAT(4X,10H(FtC!J/MIN).17X*7H(FT-IB)#/) 

111 F0RHAT(5X»F7.3»lX,F7,l,''lX»F7,2»2X#F7 . 2# 2X#?5 » 4 » 2X * F6.2) 

END 
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DESIGN PROGRAM OUTPUT SHEETS 
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LIQ HYDROGEN 

KINEMATIC V I S C 0 S I T v = 0.0 0 02 INSQ/SEC 

DENSITY = 4,3600 LB/FTCU 


UO * 

0.02 

VO = 

1,10 

RE = 

4.00 




RI = 

0 . 3 n 0 

ETA = 

79 , 500 

PT = 

1 . 525G00 

T = ') 

. 152000 


NOZZUE 

angle 

= 1.04 DEG. 






D 5 0,0150 

IN 







PT8AR 

* 17 

.32 LBP/INSQ 

HEAD 

= 4 0. 

F T 



RO 

RI 

L 

H 

N 

TOR 

hR 

VOL Fl 


(IN) 

(IN) 

(IN) 

(IN) 


( FT-t.B ) 


(FTCU/MI j) 


1 . no 

0,30 

1.00 

. 0 0 0 9 

63. 

n . o i 

0-02 

0.313 

12592. 

3.00 

0 .9Q 

3.0Q 

« 0 0 15 

182 • 

0.35 

:1 . 2 8 

4 . 688 

4197, 

5.0 0 

1.50 

5 .00 

.0019 

295. 

2 , J 5 

: 1 . 9 a 

1 6 ; 3 a 3 

2518 . 

7,00 

2.10 

7.00 

. 0 U 2 3 

405. 

6 , 51 

2*23 

371212 

1799, 

9.00 

2.7o 

9.00 

.0026 

512. 

15, 42 

4.11 

68.538 

1399 , 

11.00 

3 . 3 0 

11.00 

,0029 

6l6 . 

3J.C.5 

6 • 68 

1 1 1 4 6 0 

1149. 

13.00 

3,9c 

13.00 

, 0031 

7 i 8 - 

5 4 , 24 

10 * OU 

166 . 9q9 

969. 

15,00 

4.50 

15,00 

, 0Q33 

Sl8. 

38,33 

14-13 

235:693 

839. 

17,0 0 

5 .10 

17.00 

*0036 

$l6 . 

135.37 

19, q9 

318 .547 

741. 

19,00 

5.7q 

19.00 

.0 038 

* 1 3 . 

1 9 7 , o 3 

24.94 

4l6]lQ5 

663. 

21 ,00 

6 .30 

21.00 

. 004 0 

• 08 . 

277.37 

3i.7) 

528:963 

600 . 

23 ,00 

6,90 

23.00 

.0041 

*02 . 

373.11 

39 .42 

657 , 664 

547 . 



multiple disk turbine design study 

BASED ON 60YACK RESULTS 247 


LIQ HYDROGEN 


KINEMATIC VISCOSITY 

s 0.0002 

INSQ/SEC 




DENSITY = 

4.3600 

L3/FTCU 






UO = 

0 . 02 

vo = 

1,10 

RE 

= 4,00 




R I = 

O.3Q0 

ETA = 

79,500 

PT = 

1 . 525000 

T = 

,152000 


N02ZlE 

ANGLE 

= 1 , 04 DFG, 






D 8 0 

,0150 

IN 







PTBAR 

8 34.63 LBF/INSQ 

HEAD 

= 80 . 

FT 



RO 

Rl 

L 

H 

N 

T 0 R 

H? 

V CL FL 

RPM 

(IN) 

(IN) 

( IN ) 

(IN) 


( FT-LB) 


(FTCU/MIN) 


1.0 0 

0 .30 

1.00 

• 0 007 

64. 

0 .01 

0 • O 4 

0.375 

1 78q8 , 

3.00 

0.9q 

3.00 

•0013 

185- 

0 , 6 0 

Q • 68 

5 *. 6 5 7 

5936, 

5,00 

1.50 

5 .00 

.0016 

301. 

3,51 

2-38 

19 , 839 

3562, 

7.00 

2.10 

7.0Q 

• 0019 

414 . 

11 . 1.9 

5.42 

45 .202 

2544 , 

9,00 

2.7 0 

9.00 

.0022 

524 . 

26,55 

10-00 

6 3 . 459 

1979. 

11.00 

3 , 3 0 

11,00 

, 0 024 

632. 

52,89 

16.3 0 

136 : 013 

1619, 

13,00 

3.90 

13. 00 

, 0 026 

738. 

93,78 

24 . 46 

2 0 4 ; 0 6 3 

1370 , 

15.00 

4 . 5 il 

15,00 

• 0 029 

842 . 

1 5 3 , j 6 

34.68 

288.655 

1187, 

17,00 

5.10 

17.00 

.0030 

945. 

234,8*1 

46.83 

3 9 o 727 

1048. 

19 , no 

5 . 7 0 

19.00 

*0032 

*46 . 

3 4 3,31 

oi . 2 7 

511 ! 126 

937. 

21 , 0 0 

6,30 

21,00 

.0033 

*46 . 

4 83,01 
658,54 

77.99 

6 5 C 6 2 7 

848. 

23,0 0 

6 , 9n 

23.00 

• 0035 

*45 . 

97 • (j* 

809:943 

774, 
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LIQ HYDROGEN 

KINEMATIC VISCOSITV a 0.0002 INSO/SEC 


DENSITY S 

4,3600 

tJVFTCU 






uo = 

0 .02 

VO = 

1 .10 

RE = 

4 .00 




R I = 

0 . 3 ij o 

ETA = 

79,500 

PT = 

1 .525000 

T - r! 

. 152000 


N 0 2 Z l_ E 

ANGLE 

= 1.0 

A deg. 






D = 0 

,0150 

I M 







PTBAR 

= 51.95 lBF/INSQ 

HEAD 

= 120 . 

FT 



RO 

PI 

L 

H 

N 

TOR 

H p 

VOL fL 

RPM 

< IN > 

<IN) 

f IN) 

(IN> 


< FT-Lb ) 


( FTCU/MI ']) 


1 . 0 J 

0 . 30 

1-00 

. 0 0 0 7 

64. 

0.32 

O’O'' 

0 4 1 7 

21811. 

3.00 

0.90 

3 . 0 Q 

.0011 

186. 

0,62 

l’l-i 

6; 33 / 

7270 , 

5.00 

1 . 50 

5 ,00 

,0015 

304. 

4 ,B0 

3 • 9 g 

221164 

4362 , 

7,00 

2.10 

7,0 0 

,0017 

418. 

15,33 

9.09 

50,576 
9 3 5 0 4 

3116. 

9.00 

2,70 

9 , OC 

. 0 020 

531. 

36 , 43 

16.81 

2423 . 

ll.no 

3.30 

11.00 

, 0 022 

641. 

72,66 

27.43 

152.555 

1983. 

13.00 

3.90 

13,00 

. 0 024 

749. 

120.95 

4 l*l9 

229:u 0 

1678 . 

15,<JC 

4 . 50 

15.00 

-0 025 

855. 

210,66 

323,45 

58.32 

324 . 379 

1454 . 

1 7 , 1 o 

5.10 

17.00 

,0 027 

960. 

79-01 

439 j 449 

1283. 

l9 , no 

5 . 7 o 

19.0 0 

• 0 0 2 9 

#64 • 

473,26 

103*44 

5 7 5 3 j 4 

1148 . 

21 .00 

6.30 

21.00 

, 0 030 

#67 . 

666,31 

l3l . 76 

732 ! 845 

1039. 

23 .00 

6 . 90 

23, OC 

.0031 

#68 . 

9 )9, (i6 

164-14 

912 :9j7 

948 . 
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LIQ HYDROGEN 


kinematic viscosity 

= 0 , 

0002 

INS vj / SEC 





DENSITY = 

4,3600 

lb/ptcu 







UO = 

0 . 02 

VC = 

1.10 

RE 

= 4.00 





HI = 

n , 3 o 0 

ETA . 

79.500 

PT = 

1.525000 

T - 

u 

. 152000 


nozzle angle 

= 1.04 DEG. 







O a 

0 ,0150 

IN 








PTRAR 

s 69.26 LRF/INSQ 

head 

= 1 -) 0 . 

FT 




RO 

HI 

L 

H 

N 

TOR 

H p 


V C L . FL 

RpM 

(IN) 

(IN) 

(IN) 

(IN) 

(FT-LB) 



(FTCU/MIn) 


1.00 

0.30 

1.00 

. 0006 

64 . 

0.02 

0 

• 11 

0 .449 

25185. 

3.00 

0 ,9Q 

3 » Q 0 

.0011 

187. 

1 .02 

4 

<w 

.63 

6,811 

8395 . 

5.00 

1.50 

5.00 

.0014 

306. 

5.99 

5 

. 75 

23^965 

5037. 

7.00 

2.10 

7.00 

,0016 

421. 

19.16 

13 

•12 

54 . 74i 

3598. 

9 . 0 n 

2,70 

9.00 

.0018 

535. 

45,58 

24 

.26 

1 Q 1 ! 2 9 0 

2798 . 

11.0 0 

3 .30 

11.00 

.0020 

646 . 

90.95 

39 

.65 

165 i 383 

2290 . 

13.00 

3 . 9n 

13.00 

. 0022 

756. 

161.53 

59.58 

248'54i 

1937. 

15.00 

4.50 

15.00 

•0024 

864 . 

2 6 4, r, 4 

84.4] 

352^103 

1679. 

17.0 0 

5 .in 

l? . no 

.0025 

9?i. 

405,63 

114 

. 42 

477* 273 

1481 . 

19. 0 0 

5 . '7 o 

l 9 • 00 

• 0027 

*76. 

5 9 3,81 

: 49 

.86 

625 . i4g 

1326. 

21.0 0 

6 . 3 Q 

21.00 

• 0028 

*8 0 . 

S36 , 45 

l 9 i 

• oo 

796:713 

1199. 

23.0 0 

6.90 

23.00 

. 0029 

*83 . 

1141.71 

238 

♦ 0 <3 

9921914 

1095, 
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LI3 HYDROGEN 

KINEMATIC VISCOSITY * 0,0002 

DENSITY = 4,3600 tB/FTCU 


UO 5 

0.02 

VO = 

1,10 

RE 

RI * 

0 . 3 u 0 

ETA = 

79,500 

PT s 

NOZZLE ANqLE 

* 1.04 DEG. 


D * 

PTRAR 

0,0150 IN 
* 86.56 L3F/INSQ 

HEAD 

RO 

RI 

L 

H 

N 

<I'I> 

(IN) 

(IN) 

( IN ) 


1.00 

1 . 30 

1. OC 

. 0 006 

64 . 

3,00 

■J . 9q 

3.00 

,0010 

188. 

3.00 

1.50 

5, OC 

.0013 

3 07. 

7.00 

2.10 

7 . CO 

.0015 

424 . 

9, no 

2. 7 0 

9.00 

.0017 

538. 

11.00 

3.30 

11.00 

.0019 

65 o . 

13,00 

3.9Q 

13,0 0 

, 0 021 

761 . 

15,0 0 

4 . 5 Q 

1 5 . 0 C 

.0022 

6 7']. 

17.00 

5-10 

17.00 

,0024 

978. 

l9.nr 

5.7 0 

19.00 

. 0 025 

*85. 

21 . on 

6.30 

21.00 

. 0026 

*9n . 

23, 0 0 

6.9Q 

23. OC 

. 0028 

*95. 


i/ScC 


4,00 

1.525000 T = o. 152000 


2 C 0 . FT 


TOR 

HP 

VOL Fl 

RPM 

(CT-LB) 


(FTCJ/MI j) 


*1.03 

0 * I 4 

o'. 4 76 

28157. 

l.?l 

2-17 

7 I 227 

9386 , 

7.11 

7.6 3 

25^455 

5631. 

22.77 

17.44 

5 8 1 188 

4022. 

54,20 

32-28 

107'. 737 

3129. 

1 J 8 , 1 2 

52*74 

T 76 . 0l)6 

2560 . 

192, ?9 

79.30 

264:635 

2166 . 

3 1 4 , A 7 

112-39 

3 7 5 ! 0 7 1 

1877, 

4 63,29 

1 5 2 • 4 1 

5 0 8 6 1 5 

1656 . 

7 u 7 , 7 7 

1 9 9 ■ 7 i 

666.4-18 

1482. 

997,32 

25 4 • 61 

8491659 

1341 . 

1361.76 

317 . 4B 

*59'. 254 

1224 . 
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LU SODIUM 

KINEMATIC VISCOSITY = 0.0004 

DENSITY = 51.2000 LH/FTCU 


■ uo = 

0.02 

VO -- 

1 .10 

RE 

M I = 

0 ,3U0 

ETA r. 

79.500 

FT = 

NOZZLr. 

angle 

= 1 . L* 

4 DEG • 


D s 1 

, 0150 

IN 



PTS AR 

= 17 

.32 U-R/lNSQ 

HEAD 

Pin 

«I 

L 

H 

N 

(IN) 

(IN) 

( I'-i) 

(IN) 


1 . 10 

o.30 

1 . 0 

.0020 

59. 

3 . o J 

• .90 

3 . n.: 

. 0 0 35 

16?. 

5 . DO 

1.50 

5.0: 

. 0 045 

25 7. 

7 .00 

2.10 

7 . 0 L 

. Q () 5 3 

345 . 

9 , H. 

2 . 7 n 

9 . r o 

. 0060 

42b . 

11 • DC 

3 . 3 o 

1 1 * 0 " 

. o p 6 7 

5 o a . 

13.00 

3 . 9n 

13.0.' 

.0079 

585. 

15 • M J 

4 . 5 ;j 

1 5 • 0 Li 

• 00^8 

659. 

17.00 

5.10 

17.00 

. 0 0 8 3 

730. 

19 , DO 

5 . U) 

19.00 

. 0088 

8 o 0 . 

21.0 n 

6.30 

21.01 

.0092 

868 . 

23.00 

6 . 9 o 

2 3.0 0 

• 0 096 

934 . 


ns . / ' > „ 

H . J " 

l.S2S:'U0 T = .-.152000 


4 0. ft 


r • : R 

H P 

w OL fl 

P P y 

( F T-lB ) 


( FT CiJ/M I j ) 


0 . 0 2 

0*01 

0 .198 

36 75 . 

0.73 

0 • 1 ; 

:: . 345 

1 ? 2 6 , 

0.14 

1 , -58 

9 . 6 7 x 

735 . 


t • ? v 

21 .620 

525 . 

:< r. . 5 

2 . 3 * 

35.773 

4 03 . 

50,7 

3.74 

62:473 

334 . 

.02. 9 

5*5'. 

9?. 353 

283 . 

1 65. 76 

7 . 7 0 

129. Qj2 

245 . 

2 51 . 33 

10.3 4 

1 7 ? ; 5 9 1 

216. 

3 63,55 

13*3 • 

223:375 

193 . 

5 0 6.48 

16.87 

2 8 1 '. 3 6 5 

175. 

664.16 

20*82 

3 4 7 ; 3 1 q 

160 . 
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MULTIPLE JISK TURH I jf DESIGN SH'l.* 
BASED 3N BOYACK RESULTS 


LIQ sodium 

KINEMATIC VISCOSIT v = 0.0004 iNS./bcr 

DENSITY = 51,2000 lB/PTCU 


uo = 

0.02 

VO = 

1 ,10 

RE = 

< . c w 




RI = 

O . 3 () o 

ETA = 

7 o , 5 0 0 

PT = 

i.'it 0 n 

T - ii, 

,152033 


VOZZLE 

ANGLE 

= J.i- 

4 DE G • 






n = j 

,0150 

IN 







P T « A ^ 

= 34,65 LBF/INSQ 

heap 

= * o , 

PT 



RD 


L 

H 

IN 

tqR 

H P 

VCL V L 

R P *>' 

( I N ) 

(IN) 

( I "J ) 

(IN) 


CFT-LE) 

( 

( P T C U / M I \j ) 


l . no 

O .30 

l.O? 

.3317 

6 0 . 

.j , 3 

"•02 

r .24 j 

5197 . 

■5 , o 0 

0.90 

3.00 

.0 029 

167 . 

1 » 2 7 

" . 4 c 

3,487 

1732 . 

5 . n (! 

1.50 

5.30 

. 0 3 3 8 

266 . 

7.23 

1.42 

11.939 

1 0 39 . 

7 . U- 

2 . 10 

7 . 3 0 

, 0 n 4 5 

3 6() . 

22, -,4 

3.20 

26 , 7q3 

7 4 2 , 

9 , i 0 

2 , 7o 

9 » to 

. 0 3 5 1 

449 . 

> 2 , • 4 

5 .82 

45,563 

57 7 . 

11 . 't ■ 

3 . 3(j 

1 1 . U 

. 0 0 5 6 

534 . 

1)4, ;> 

9 . 36 

7S : 116 

4 7 2. 

13 , 

3 . 9 o 

1 3 . I 0 

• 0 361 

6i6 . 

1 o 2 .47 

13.8 9 

115.864 

4 0 0 . 

IS , '10 

4.50 

15 . 3 j 

• 0365 

696 . 

994,73 

19-44 

162.224 

246 , 

17.00 

5.10 

17.3 j 

. 0 3 7 0 

774 . 

4 4 8, 4 

26-08 

21 7 i 553 

306 . 

19,00 

5 . /C 
5.30 

19.0 0 

• 0 3 74 

85 0 . 

'■4 9,3 0 

53.8 - 

232; 185 
3561338 

274 , 

21.00 

21.30 

.0 377 

924 . 

5 1 6 , i) 4 

4 2-74 

247. 

2 3.03 

6.90 

23.00 

.0 3 81 

996 . 

12 2 7.33 

5 2.79 

4 4 3 ’.423 

226 . 
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LIQ SODIUM 

KINEMATIC VISCCSIT' = 0.0004 iMSwSZr: 

DENSITY = 51,20 0 0 i.b/ETCU 


uo = 

o . 02 

V0 - 

1.10 

RE = 

4 ■ J 1 




RI = 

o .3:0 

eta -- 

79.500 

PT = 

1 . 523. 0i 

T : J 

. 152Q00 


Li I 
-J 
•Vi 
TV, 

o 

A i'! G L E 

= 1 . 

a JF:G . 






U = o 

. 0150 

IM 







PTBAR 

= 51 

• 95 I NSQ 

HEAD 

- 12 0. 

R T 



RO 


L 

H 

!•; 

1 C R 

H p 

'■cl r l 

R P M 

( I 'i ) 

( I N ) 

(IN) 

(IN) 


( v T-l.L' ) 


( RTLU/M I 0 ) 


1 . OC 

■ : . 3 0 

l . n o 

. 0 U5 

61. 

0.1 4 

C • ? 5 

0 .269 

6365 . 

3 . n c 

0.90 

3 . C 

.0026 

170. 

1.75 

;■ . 7 1 . 

3.921 

2122 . 

5 . 0 C 

1.50 

5 . n 

. o '< 3 4 

272. 

i c . 0 0 

2-47 

13.474 

1275. 

7 . ;1 n 

2.10 

7. [)p 

. 0 0 4 o 

368 . 

31.79 

5 . 4 J 

30.221 

9 0 9 . 

9.-0 

2. /:) 

9 . 0 L 

. 0 ' 1 4 6 

460 . 

7 3. V 

9.90 

55.082 

7 0 7 . 

11.0 0 

3.30 

11 • Ul 

• 0 ‘1 5 1 

548 . 

j 44 .; 9 

15.9' 

*6 . 776 

579 . 

13* 00 

3 • 9 o 

13.00 

. 0 !;55 

634 . 

. * 4 , 

2 3 • 7 , 

j 51.895 

4 9 r , 

15 . '1 

4 . 5 0 

15.0 0 

.0059 

^ 1 7 . 

4U.I9 

3 8*2 p 

104.944 

424 . 

17.0 0, 

5.10 

17.00 

. 0 lift 3 

798. 

'• 2 6 . -*2 

44.65 

2 4 8 . 3 8 5 

374 . 

19 . P 0 

5 . / 0 

19.0 0 

. 0 0 6 7 

078. 

' ■ 9 , 1 9 

57.99 

322.524 

335, 

21.'!'" 

ft . 30 

2 1 . j 3 

. 0 i' 7 0 

9 55. 

1 , 7 0 

73 . 37 

407.787 

3 0 3 . 

23.0 0 

6.90 

23 .00 

. 0 0 7 3 

*31 . 

J 7 2 1 , *i a 

9. .70 

5 0 4 . 455 

27 7. 
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LID SODIUM 

XIMEUTIC V I 5 C 0 3 I T v = 3,0004 153 j / 3L J 

DENSITY = 51,200 0 u'VETCU 


00 = 

0.02 

VO = i 

.10 

RE = 

4 . 0 1 




SI = 

i . 3 j n 

ETA -.71 

.500 

PT = 

1 . 5 .2 5 i 0 

T = . 

i 1 5 2 0 1 0 


NOZZLE 

« 

A Ur,L E 

= 1.04 

DEG. 






D - v 

, ■"* 15 0 

I N 







ptrar 

= 69. 

2 6 L B T / l 

us a 

head 

1 ' ' j . 

► T 



RO 

9 I 

L 

H 

r i 

'{ 

T y-i 

H ^ 

VOL LL 


< I’D) 

(IN) 

( I N ) 

( IN ) 


< " T-lB ) 

( 

E T U J / M I 1 ) 


1 , 0 0 

0 . 3 0 

1.00 . 

0 ! 1 4 

61. 

-i. 15 

: - a/ 

9 . 2 9 7 

7 5 4 9. 

3 . On 

J .9Q 

3.10 

0 0 25 

172. 

2 . 19 

1-02 

4 '. 258 

2 4 5 J . 

5,0 0 

1.50 

5 . 0 o 

0132 

275. 

. 2 . •> 6 

3 • 5 7 

14 .667 

14 7 0 . 

7 . 0 0 

2.10 

7 . 0 o 

3 0 38 

37 3. 

39,53 

7.9) 

32 . 95 9 

1 C 5 0 . 

9 , 0 n 

2.7o 

9.10 

0 3 43 

467. 

'> 2 . n 

14.47 

60.165 

617. 

11-00 

3 . 3 1 1 

11-00 

0 047 

558. 

1»2 , -'8 

2 5-2* 

9 7 . n 9 6 

663 . 

13.03 

3 , 9 0 

1 3 • J 0 . 

0 0 5 1 

646 . 

3 5 1 , .5 

3 4.6/ 

1 4 4 . 4 2 J 

665 . 

1 5 . n n 

4.5J 

15.00 • 

i ; 5 5 

732- 

6 2.1 . 3 

4 3 • 6 i 

2 u 2 . 7 1 1 

4 9 (j , 

17,00 

5 . 10 

17.00 . 

0 3 59 

815. 

79J , 94 

6 L . 3 , 

2 7 2 . 4 c 6 

432. 

19.1) 0 

3 . / n 

19.00 

0762 

397. 

1152.71 

5 4.8 7 

354,124 

33/. 

21 . no 

6 , 3 0 

21.00 

0 r,6 5 

976. 

1 6 1 ? , i. 3 

1 0 7 , 4 / 

448.032 

35 0 . 

23.00 

6.9c 

23.00 

1068 

#55 . 

2115 . 7j. 

132.98 

5 5 4 ; 6 9 .4 

32 0 . 
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LID SODIUM 

KINEMATIC VISCOSITY = 0.0304 I NSU/ S.;U 

DEN S I r Y = 51.2000 LD/TTCU 


uo = 

0 . C 2 

VO = 1 

.10 

RE = 

4 . 0 1 




RI = 

0 . 3 u 0 

ETA =79 

.500 

PT = 

1.525 0 0 1 

T = 0 

. 1520 0 0 


'-■07ZLE 

ANGLE 

= 1 . I 4 

') E G . 






ns o 

. 0150 

IN 







PTRA.R 

= 86 

.56 let/i 

r ;SQ 

HEAD 

= 2l J 

FT 



Rn 

Si 

L 

H 

K- 

T J R 

H p 

VOL r|_ 

R P M 

(IN) 

( r. ) 

< 1 N ) 

(IN) 


( r I ■" L. B 


( r T c U / M i - : ! 


1 . 0 0 

o . 3 0 

1 . n 0 

no 13 

61. 

0 . ;6 

- . P 'v 

0 . 3 j f 

621 7 . 

3 . no 

J . 9 0 

3 . U 0 

n j23 

173. 

2.61 

u-3'' 

4 ; 5 3 7 

2739 , 

5 . on 

1.50 

5 . OU • 

0 0 3 n 

278. 

\ 4 . u 9 

4.69 

1 5 6 j i 

1643 . 

7 . 0 'i 

2.10 

7 , no 

0 , 1 3 6 

377. 

4 7.24 

x ). r >6 

33 . 2 33 

1174 . 

9.00 

2 .70 

9.00 

n 0 4 0 

473. 

mi; \ , i 

19 . 3 ■' 

6 4.39:} 

9.13 . 

11.00 

6 . 3n 

1 1 . 0 C 

n ■ 1 4 5 

565. 

v) y g 4 k 

3 . 1 / 

104 . ni5 

7 47 . 

13.00 

3.9;j 

13-00 . 

0 ,0 4 8 

655. 

365 7 

4 6-4 0 

154,844 

632. 

1 5 . iin 

4.50 

15*00 

0 ■ j 5 2 

743. 

0 2 l . - 2 

’n • 1" 

217:5]? 

546 . 

17.00 

5. 10 

1 7 . ) U . 

0 0 55 

828. 

952 . ><2 

6 7.67 

2 9 p . 8 5 t 

463 . 

1 9 . 0 ) 

5.7 0 

19. oc 

^ n59 

911. 

1358..- I 

1 1 4 • 0 1 

<60 . 4 

432. 

21.00 

0.33 

21. UP 

i 1 0 6 2 

993. 

1 9 3 7 , 5 / 

144 • 31 

481.67/ 

391. 

23,0! 

6 .9^ 

23.0 0 . 

0 0 6 4 

*73. 

■'6 26 . 3 

1 78.7V 

5 9 6.6 3 n 

357. 



y UL T l PLfc DISK TURiUNF 'ESIGN STUO v 
BASED ON bHYACK RESULT.-' 


GLYCERIN 


KINEMATIC VISCOSITY 

= 0 . 

0323 

I Ns N/ SEC 




DENSIT 

v = 79 

.3000 

LRATCU 






UO = 

0 . ') c 

VO = 

1 . 1 C 

PE 

= 4,00 




RI = 

o . 3 ' . 0 

ETA = 

79.500 

PT = 

1 .5280 p ii 

T = 

, 1520PC 


NOZZlF. 

ANGLE 

s 1 . 0 

4 DEG. 






D = 0 

. D15C 

IN 







PTPAR 

= 17. 

32 LET/IrjSQ 

HEAD 

= . 

P T 



RO 

RI 

L 

H 

N 

TOP 

pR 

w D L r L 

RpM 

UN> 

( I s ) 

( IM 

(IN) 

f FT-Lb) 


( FT CU/ M I ' > 


1 . n 0 

f, .30 

i • 0 t* 

• 0 f- 0 4 

28. 

0 ■ • 1 8 

(■ > 0 5 

C . 7 7 2 . 

2953 . 

3. CD 

t.90 

3.00 

. 0 354 

60 . 

2 . 7 3 

C • 5 1 

3.529 

984 . 

5 , or, 

1.50 

5.00 

.0457 

82 . 

a 3 . 1 4 

1 > 5 2 

25.395 

591. 

7,00 

2 . 1C 

7. DC 

. 0 5 4 1 

101 • 

38.02 

5 • 1 •' 

51.754 

422 . 

9 , fj n 

2.70 

9 . ') o 

.05J.3 

118 . 

54 . ; 4 

5 • 

8 7.804 

3 2 5. 

11-5 0 

3.3 0 

11*00 

. p P 7 8 

133 . 

1.85 . 05 

8.01 

1 3 3 . 6 7 2 

266 . 

13.00 

3.90 

1 3 . 0 C 

. p 7 3 7 

147. 

Z 82 , 55 

11-35 

169.449 

2 2 7. 

15,0 0 

4.50 

15.0 0 

.0 792 

159. 

4 l- t . n 

15.3 0 

255.206 

197. 

17,0 i : 

5.10 

17.0Q 

. 0 4 3 

171. 

39 9,1-6 

1 9 . d 4 

330.996 

174. 

19,00 

5.70 

19. BO 

. 0 ' 9 1 

183. 

844 . 36 

24.9 S 

4 1 6 . 8 6 6 

155 . 

21 . nc 

6 .30 

21.00 

. 0 > 3 7 

193. 

.1148.13 

3(1.74 

5 1 2 i 8 3 2 

i 4 1 . 

23 . OC 

6 . 9 o 

23 . OC 

.0980 

203 . 

1517, 72 

3 7.10 

618,98 9 

128. 
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GLYCERIN 


KINEMATIC V I s C 0 S I T y 

= 0 , 

0323 

I N$0/S£L 




dens I ty 

= 79.3000 

lb/ctcu 






UO 5 0 

. 02 

VO = 

1 . 1 ” 

RE 

- ‘i . 0 w 


192000 


PI « n 

. 3 on 

ETA = 

79,3 .1 0 

FT = 

i. .*25 if..' 

T * ’ , 


NOZZLE 

angle 

= 1 , 0 

4 DEG* 






D 3 ' . 

H5n 

IN 







PJHAH = 

3 4 

.93 LBF/I jSU 

READ 

^ 1 ■' f 

T 



RO 

RI 

L 

H 

N 

T i , 4 R 

n p 

VOL |-L 

9P‘4 

< ri > 

< I >' > 

I I N ) 

< IN 

) 

< FI -LB > 

< 

R T 0 J / 4 1 


1.0 0 

r ' . 3 0 

1 . O'.i 

* 0 1 7 2 

31. 

0 * 

0 * 1 2 

1 .>U 9 

4 i ’’6 . 

3.00 

0 . V o 

3,00 

• [? V ° 8 

67. 

9 , l 6 

1*37 

11 . 4 19 

13 92. 

5 ,00 

1.90 

9.0C 

* 0 < R 4 

94 . 

? 5 . 8 6 

4.11 

34.3 M 

6.35 . 

7 , OR 

2*10 

7 . 0 0 

• 0 *99 

116 . 

74, ,0 

1 . 4 3 

7 0.3 0 3 

5 97 . 

9,00 

2 * 7 n 

9*00 

• ()9l6 

135 . 

..62,4 5 

j 4.39 

1 1 9 . 722 

4 6 4. 

11-00 

3 * 3 n 

11*0 0 
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GLYCERIN 

KINEMATIC VISCOSITY = 0.0323 INgC/SEC 

DENSITY = 79,3000 LR/FTCU 
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DENSITY = 4,3600 LB/TTCU 

UO 5 0,02 VO = 1.10 RE = 

PI s 0,300 ETA =79,500 PT = 


NOZZLE 

ANGLE 

= 1.04 

DEG . 

D = 0 

,0150 

IN 


PTPAR 

= 86 

.58 LET/InSO HEAD 

RO 

R 1 

L 

h N 

(IN) 

(IN) 

(IN) 

( IN ) 

2.00 

0.60 

2. OC 

. o o o a 12 6. 


OFF DESIGN OPERATION 
WITH HEAD VARYING 


VOL FL 

RpM 

hP 

TOR 

(FTCU/MIN ) 


(ft-lb > 

2.906 

15422.4 

1.05 

0 . 36 

2.845 

15097.7 

G . 98 

0.34 

2.782 

14765.9 

0.92 

0.33 

2.719 

14426.4 

0 . 86 

0.31 

2.653 

14078.7 

0.80 

0.30 

2.586 

13722.2 

0.74 

0 . ?8 

2.517 

13356.2 

0.68 

0.27 

2.446 

12979,9 

0.62 

0 . 25 

2.3 73 

12592.4 

fi .57 

0 . 24 

2.298 

12192 . 5 

C.52 

0 . 22 

2.220 

11779.1 

0.47 

0.21 


4,00 

1,525000 T = ',152000 


2 0 0. FT 


TOR 

1 i P v 0 L F L PPM 

FT-LB ) 

( F T r U / M I N ) 

0 .30 

0*80 2.65 i)079 


F.Ta 

T p A f- 

793? 

1.'3 , 9f, 

7932 

99,57 

7932 

^ 4 

793? 

9 0.91 

7932 

36 , 36 

793? 

5 2. 25 

793? 

77.9? 

793? 

7 3,59 

7932 

6 9,2ft 

7932 

64,94 

793? 

H ' . 6 i 



MULTIPLE DISK TURBINE DESIGN STUDY 
BASED ON BOYACK RESULTS 


262 


LIQ SODIUM 
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GLYCERIN 

KINEMATIC VISCOSITY « 0.0323 INSfJ/SEC 
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